Elliptical polynomial-based message authentication code

ABSTRACT

The elliptic-polynomial based Message Authentication Code (MAC) provides MAC generation methods based on the elliptic polynomial discrete logarithm problem. It is well known that an elliptic polynomial discrete logarithm problem is a computationally “difficult” or “hard” problem. The methods use both an elliptic polynomial and its twist, even if the polynomial and its twist are not isomorphic. Since both the polynomial and its twist are used, multiple x- and y-coordinates can be used to embed bit strings into a point that satisfies the elliptic polynomial, and the embedding process is non-iterative, so that the time required to embed the bit string is independent of the bit string content.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to computerized cryptographic methods for communications in a computer network or electronic communications system, and particularly to an elliptic polynomial-based message authentication code.

2. Description of the Related Art

In recent years, the Internet community has experienced explosive and exponential growth. Given the vast and increasing magnitude of this community, both in terms of the number of individual users and web sites, and the sharply reduced costs associated with electronically communicating information, such as e-mail messages and electronic files, between one user and another, as well as between any individual client computer and a web server, electronic communication, rather than more traditional postal mail, is rapidly becoming a medium of choice for communicating information. The Internet, however, is a publicly accessible network, and is thus not secure. The Internet has been, and increasingly continues to be, a target of a wide variety of attacks from various individuals and organizations intent on eavesdropping, intercepting and/or otherwise compromising or even corrupting message traffic flowing on the Internet, or further illicitly penetrating sites connected to the Internet.

Encryption by itself provides no guarantee that an enciphered message cannot or has not been compromised during transmission or storage by a third party. Encryption does not assure integrity due to the fact that an encrypted message could be intercepted and changed, even though it may be, in any instance, practically impossible, to cryptanalyze. In this regard, the third party could intercept, or otherwise improperly access, a ciphertext message, then substitute a predefined illicit ciphertext block(s), which that party, or someone else acting in concert with that party, has specifically devised for a corresponding block(s) in the message. The intruding party could thereafter transmit the resulting message with the substituted ciphertext block(s) to the destination, all without the knowledge of the eventual recipient of the message.

The field of detecting altered communication: is not confined to Internet messages. With the burgeoning use of stand-alone personal computers, individuals or businesses often store confidential information within the computer, with a desire to safeguard that information from illicit access and alteration by third parties. Password controlled access, which is commonly used to restrict access to a given computer and/or a specific file stored thereon, provides a certain, but rather rudimentary, form of file protection. Once password protection is circumvented, a third party can access a stored file and then change it, with the owner of the file then being completely oblivious to any such change.

In order to ensure message integrity, systems that depend on a shared secret key often rely upon message authentication codes (MACs). The MAC is uses the secret key and the message text to produce a bit string unique to the key and the text of the message, the bit string being appended to the message. The MAC bit string is usually compressed to a fixed number of bits, so that it is much shorter than the message text. When a recipient is able to use his secret key and the text message to generate a MAC that is identical to the MAC appended to the message, the recipient can be confident in the integrity of the message, i.e., that the message has not been altered.

A number of encryption techniques have been used to generate MACs. Some MACs are generated using block ciphers, such as the MAC made with DES. Such MACs, however, are vulnerable to brute force attacks. Other MACs are generated using hash functions, such as MD5, RIPEMD-160, SHA-1, etc. One popular MAC is HMAC, which combines a secret key with a non-keyed hash function, such as SHA-1 or MD5. However, such MACs rely upon the invulnerability of the hash function for their security against various forms of attack. However, the security of such hash functions against attacks has come into question because of successful attacks and evidence of their vulnerability to collisions. None of the current MAC generation algorithms incorporate a hash function that is based upon a mathematically hard problem, and particularly not the elliptic polynomial discrete logarithm problem.

Thus, an elliptic polynomial-based message authentication code solving the aforementioned problems is desired.

SUMMARY OF THE INVENTION

The elliptic-polynomial based Message Authentication Code (MAC) provides MAC generation methods based on the elliptic polynomial discrete logarithm problem. It is well known that an elliptic polynomial discrete logarithm problem is a computationally “difficult” or “hard” problem. The methods use both an elliptic polynomial and its twist, even if the polynomial and its twist are not isomorphic. Since both the polynomial and its twist are used, multiple x- and y-coordinates can be used to embed bit strings into a point that satisfies the elliptic polynomial, and the embedding process is non-iterative, so that the time required to embed the bit string is independent of the bit string content.

These and other features of the present invention will become readily apparent upon further review of the following specification and drawings.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The elliptic polynomial-based Authentication Code (MAC) provides MAC generation methods based on the elliptic polynomial discrete logarithm problem. It is well known that an elliptic polynomial discrete logarithm problem is a computationally “difficult” or “hard” problem.

The MACs to be described below use elliptic polynomials in their generation, where different elliptic polynomials are used for different blocks of the same plaintext. Particularly, the MACs use an elliptic polynomial with more than one independent x-coordinate. More specifically, a set of elliptic polynomial points are used which satisfy an elliptic polynomial equation with more than one independent x-coordinate which is defined over a finite field F having the following properties: One of the variables (the y-coordinate) has a maximum degree of two, and appears on its own in only one of the monomials; the other variables (the x-coordinates) have a maximum degree of three, and each must appear in at least one of the monomials with a degree of three; and all monomials which contain x-coordinates must have a total degree of three.

The group of points of the elliptic polynomial with the above form is defined over additions in the extended dimensional space, and, as will be described in detail below, the method makes use of elliptic polynomials where different elliptic polynomials are used for different blocks of the same plaintext.

The particular advantage of using elliptic polynomial cryptography with more than one x-coordinate is that additional x-coordinates are used to embed extra message data bits in a single elliptic point that satisfies the elliptic polynomial equation. Given that nx additional x-coordinates are used, with nx being greater than or equal to one, a resulting elliptic point has (nx+1) x-coordinates, where all coordinates are elements of the finite field F. The number of points which satisfy an elliptic polynomial equation with nx additional x-coordinates defined over F and which can be used in the corresponding cryptosystem is increased by a factor of (#F)^(nx), where # denotes the size of a field.

Through the use of this particular method, security is increased through the usage of different elliptic polynomials for different message blocks during the generation of a message authentication code. Further, each elliptic polynomial used for each message block is selected at random, preferably using an initial value and a random number generator.

Given the form of the elliptic polynomial equation described above, the elliptic polynomial and its twist are isomorphic with respect to one another. The method uses an embedding technique, to be described in greater detail below, which allows the embedding of a bit string into the x-coordinates of an elliptic polynomial point in a deterministic and non-iterative manner when the elliptic polynomial has the above described form. This embedding method overcomes the disadvantage of the time overhead of the iterative embedding methods used in existing elliptic polynomial.

The difficulty of using conventional elliptic polynomial cryptography to develop MACs typically lies in the iterative and non-deterministic method needed to embed a bit string into an elliptic polynomial point, which has the drawback of the number of iterations needed being different for different bit strings which are being embedded. As a consequence, different calculation times are required for different blocks of bit strings. Such a data-dependant generation time is not suitable for generating MACs, which require data independent encryption time. Further, with regard to iterative and non-deterministic methods in conventional elliptic polynomial, cryptography, given an elliptic polynomial defined over a finite field that needs N-bits for the representation of its elements, only ((nx+ny+1)N−L) bits of the message data bits can be embedded in any elliptic polynomial point.

The isomorphic relationship between an elliptic polynomial and its twist, which is obtained as a result of the given form of the elliptic polynomial equation, ensures that any bit string whose equivalent binary value is an element of the underlying finite field has a bijective relationship between the bit string and a point which is either on the elliptic polynomial or its twist. This bijective relationship allows for the development of the elliptic polynomial MACs to be described below.

In the conventional approach to elliptic polynomial cryptography, the security of the resulting cryptosystem relies on breaking the elliptic polynomial discrete logarithm problem, which can be summarized as: given the points k(x_(0,B), x_(1,B), . . . , x_(nx,B), y_(B)) and (x_(0,B), x_(1,B), . . . , x_(nx,B), y_(B)), find the scalar k.

Further, projective coordinates are used at the sending and receiving entities in order to eliminate inversion or division during each point addition and doubling operation of the scalar multiplication. It should be noted that all of the elliptic polynomial cryptography-based MACs disclosed herein are scalable.

In the following, with regard to elliptic polynomials, the “degree” of a variable u^(i) is simply the exponent i. A polynomial is defined as the sum of several terms, which are herein referred to as “monomials”, and the total degree of a monomial u^(i)v^(j)w^(k) is given by (i+j+k). Further, in the following, the symbol ε denotes set membership.

One form of the subject elliptic polynomial equation with more than one x-coordinate and one or more y-coordinates is defined as follows: the elliptic polynomial is a polynomial with more than two independent variables such that the maximum total degree of any monomial in the polynomial is three; at least two or more of the variables, termed the x-coordinates, have a maximum degree of three, and each must appear in at least one of the monomials with a degree of three; and at least one or more variables, termed the y-coordinates, have a maximum degree of two, and each must appear in at least one of the monomials with a degree of two.

Letting S_(nx) represents the set of numbers from 0 to nx (i.e., S_(nx)={0, . . . , nx}), and letting S_(ny) represents the set of numbers from 0 to ny (i.e., S_(ny)={0, . . . , ny}), and further setting (nx+ny)≧1, then, given a finite field, F, the following equation defined over F is one example of the polynomial described above:

$\begin{matrix} {{{{\sum\limits_{k \in S_{ny}}\;{a_{1k}y_{k}^{2}}} + {\sum\limits_{k,{l \in S_{ny}},{l \neq k}}\;{a_{2{kl}}y_{k}y_{l}}} + {\sum\limits_{k \in S_{ny}}\;{a_{3k}y_{k}}} + {\sum\limits_{k,{l \in S_{ny}},{i \in S_{nx}}}\;{c_{1{kli}}y_{k}y_{l}x_{i}}} + {\sum\limits_{{k \in S_{ny}},{l \in S_{nx}}}\;{c_{2{kl}}y_{k}x_{l}}} + {\sum\limits_{{k \in S_{ny}},l,{i \in S_{nx}}}\;{c_{3{kli}}y_{k}x_{l}x_{i}}}} = {{\sum\limits_{l \in S_{nx}}\;{b_{1l}x_{l}^{3}}} + {\sum\limits_{l,{k \in S_{nx}},{l \neq k}}\;{b_{2{lk}}x_{l}^{2}x_{k}}} + {\sum\limits_{l,{k \in S_{nx}}}\;{b_{3{lk}}x_{l}x_{k}}} + {\sum\limits_{k \in S_{nx}}\;{b_{4k}x_{k}}} + b_{c}}},} & (1) \end{matrix}$ where a_(1l), a_(2kl), a_(3k), c_(1lki), c_(2kl), c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k) & b_(c)εF.

Two possible examples of equation (1) are y₀ ²=x₀ ³+x₁ ³+x₀x₁ and y₀ ²+x₀x₁y₀+y₀=x₀ ³+x₁ ³+x₀ ²x₁+x₀x₁ ²+x₀x₁+x₀+x₁.

With regard to the use of the elliptic polynomial equation in the addition of points of an elliptic polynomial with more than one x-coordinate and one or more y-coordinates, we may examine specific coefficients a_(1k), a_(2kl), a_(3k), c_(1lki), c_(2kl), c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k) & b_(c)εF for F, wherein a set of points EC^(nx+ny+2) is defined as the (nx+ny+2)−tuple (x₀, x₁, . . . , x_(nx), y₀, y₁, . . . y_(ny)), where x_(i), Y_(k)εF, iεS_(nx) and kεS_(ny). This set of points comprises solutions of F, although excluding the point (0, 0, . . . , 0) and the point at infinity, (x_(0,I), x_(1,I), . . . , x_(nx,I), y_(0,I), y_(1,I), . . . , y_(ny,I)).

The rules for conventional elliptic polynomial point addition may be adopted to define an additive binary operation, “+”, over EC^(nx+ny+2), i.e., for all (x _(0,1) ,x _(1,1) , . . . , x _(nx,1) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))εEC ^(nx+ny+2) and (x _(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2))εEC ^(nx+ny+2), the sum: (x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=(x _(0,1) ,x _(1,1) , . . . , x _(nx,1) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))+(x _(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y _(0,2) y _(1,2) , . . . , y _(ny,2)) is also (x _(0,3) ,x _(1,3) , . . . x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))εEC ^(nx+ny+2).

As will be described in greater detail below, (EC^(nx+ny+2), +) forms a pseudo-group (p-group) over addition that satisfies the following axioms:

-   -   (i) There exists a set (x_(0,I), x_(1,I), . . . , x_(nx,I), . .         . , y_(0,I), y_(1,I), . . . , y_(ny,I))εEC^(nx+ny+2) such that         (x₀, x₁, . . . , x_(nx), y₀, y₁, . . . ,         y_(ny))+(x_(0,I)x_(1,I), . . . x_(nx,I), y_(0,I), y_(1,I), . . .         , y_(ny,I))=(x₀, x₁, . . . , x_(nx), y₀, y₁, . . . , y_(ny)) for         all (x₀, x₁, . . . , x_(nx), y₀, y₁, . . . ,         y_(ny))εEC^(nx+ny+2);     -   (ii) for every set (x₀, x₁, . . . , x_(nx), y₀, y₁, . . . ,         y_(ny))εEC^(nx+ny+2), there exists an inverse set, −(x₀, x₁, . .         . , x_(nx), y₀, y₁, . . . , y_(ny))εEC^(nx+ny+2), such that (x₀,         x₁, . . . , x_(nx), y₀, y₁, . . . , y_(ny))−(x₀, x₁, . . . ,         x_(nx), y₀, y₁, . . . , y_(ny))=(x_(0,I), x_(1,I), . . . ,         x_(nx,I), y_(0,I), y_(1,I), . . . , y_(ny,I))     -   (iii) the additive binary operation in (EC^(nx+ny+2), +) is         commutative, and the p-group (EC^(nx+ny+2), +) forms a group         over addition when:     -   (iv) the additive binary operation in (EC^(nx+ny+2), +) is         associative.

Prior to a more detailed analysis of the above axioms, the concept of point equivalence must be further developed. Mappings can be used to indicate that an elliptic point represented using (nx+1) x-coordinates and (ny+1) y-coordinates, (x₀, x₁, . . . , x_(nx), y₀, y₁, . . . , y_(ny)), is equivalent to one or more elliptic points that satisfy the same elliptic polynomial equation, including the equivalence of an elliptic point to itself.

Points that are equivalent to one another can be substituted for each other at random, or according to certain rules during point addition and/or point doubling operations. For example, the addition of two points (x_(0,1), x_(1,1), . . . , x_(nx,1), y_(0,1), y_(1,1), . . . , y_(ny,1)) and (x_(0,2), x_(1,2), . . . , x_(nx,2), y_(0,2), y_(1,2), . . . , y_(ny,2)) is given by: (x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=(x _(0,1) ,x _(1,1) , . . . , x _(nx,1) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))+(x_(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2))

If the point (x_(0,1)″, x_(1,1)″, . . . , x_(nx,1)″, y_(0,1)″, y_(1,1)″, . . . , y_(ny,1)″) is equivalent to the point (x_(0,1), x_(1,1), . . . , x_(nx,1), y_(0,1), y_(1,1), . . . , y_(ny,1)), then the former can be substituted for (x_(0,1), x_(1,1), . . . , x_(nx,1), y_(0,1), y_(1,1), . . . , y_(ny,3)) in the above equation in order to obtain: (x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=(x _(0,1) ″,x _(1,1) ″, . . . , x _(nx,1) ″,y _(0,1) ″,y _(1,1) ″, . . . , y _(ny,1))+(x _(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2))

Mappings that are used to define equivalences can be based on certain properties that exist in elliptic polynomial equations, such as symmetry between variables. As an example, we consider the point (x₀,x₁,y₀) that satisfies the equation y₀ ²=x₀ ³+x₁ ³+x₀x₁. The equivalent of this point may be defined as (x₁,x₀,−y₀).

With regard to the addition rules for (EC^(nx+ny+)2, +), the addition operation of two points (x_(0,1), x_(1,1), . . . , x_(nx,1), y_(0,1), y_(1,1), . . . , y_(ny,1))εEC^(nx+ny+2) and (x_(0,2), x_(1,2), . . . , x_(nx,2), y_(0,2), y_(1,2), . . . , y_(ny,2))εEC^(nx+ny+2), otherwise expressed as: (x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=(x _(0,1) ,x _(1,1) , . . . , x _(nx,1) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))+(x _(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2)) is calculated in the following manner. First, a straight line is drawn which passes through the two points to be added. The straight line intersects EC^(nx+ny+2) at a third point, which we denote (x_(0,3)′, x_(1,3)′, . . . , x_(nx,3)′, y_(0,3)′, y_(1,3)′, . . . , y_(ny,3)′)εEC^(nx+ny+2). The sum point is defined as (x_(0,3), x_(1,3), . . . , x_(nx,3), y_(0,3), y_(1,3), . . . , y_(ny,3))=−(x_(0,3)′, x_(1,3)′, . . . , x_(nx,3)′, y_(0,3)′, y_(1,3)′, . . . , y_(ny,3)′).

From the above definition of the addition rule, addition over EC^(nx+ny+2) is commutative, that is: (x _(0,1) ,x _(1,1) , . . . , x _(nx,1) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))+(x _(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2))=(x _(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2))+(x _(0,1) ,x _(1,1) , . . . , x _(nx,1) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1)) for all (x_(0,1), x_(1,1), . . . , x_(nx,1), y_(0,1), y_(1,1), . . . , y_(ny,1))εEC^(nx+ny+2) and for all (x_(0,2), x_(1,2), . . . , x_(nx,2), y_(0,2), y_(1,2), . . . , y_(ny,2))εEC^(nx+ny+2). This commutivity satisfies axiom (iii) above.

There are two primary cases that need to be considered for the computation of point addition for (EC^(nx+ny+2), +): (A) for at least one jεS_(nx), x_(j,1)≠x_(j,2); and (B) for all jεS_(nx),x_(j,1)=x_(j,2)=x_(j,0). Case B includes three sub-cases:

-   -   i. for all kεS_(ny) y_(k,1)=Y_(k,2), that is:         (x _(0,1) ,x _(1,1) , . . . , x _(nx,1) ,y _(0,1) ,y _(1,1) , .         . . , y _(ny,1))=(x _(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y         _(0,2) ,y _(1,2) , . . . , y _(ny,2)),     -   which corresponds to point doubling;     -   ii. for kεS_(ny) & k≠0, y_(k,1)=y_(k,2), and where y_(0,1) &         y_(0,2) are the roots of the following quadratic equation in y₀:

${{{a_{10}y_{0}^{2}} + {\sum\limits_{{k \in S_{ny}},{k \neq 0}}\;{a_{1k}y_{k,1}^{2}}} + {y_{0}\left\{ {{\sum\limits_{{k \in S_{ny}},{k \neq 0}}\;{a_{2k\; 0}y_{k,1}}} + {\sum\limits_{{l \in S_{ny}},{l \neq 0}}\;{a_{20l}y_{l,1}}}} \right\}} + {\sum\limits_{k,{l \in S_{ny}},{l \neq k},{{{l\&}k} \neq 0}}\;{a_{2{kl}}y_{k,1}y_{l,1}}} + {a_{30}y_{0}} + {\sum\limits_{{k \in S_{ny}},{k \neq 0}}\;{a_{3k}y_{k,1}}} + {y_{0}^{2}{\sum\limits_{i \in S_{nx}}\;{c_{100i}x_{i,1}}}} + {y_{0}\left\{ {{\sum\limits_{{k \in S_{ny}},{i \in S_{nx}}}\;{c_{1k\; 0i}y_{k,1}x_{i,1}}} + {\sum\limits_{{l \in S_{ny}},{i \in S_{nx}}}\;{c_{10{li}}y_{l,1}x_{i,1}}}} \right\}} + {\sum\limits_{k,{l \in S_{ny}},{{{l\&}k} \neq 0},{i \in S_{nx}}}\;{c_{1{kli}}y_{k,1}y_{l,1}x_{i,1}}} + {y_{0}{\sum\limits_{l \in S_{nx}}\;{c_{20\; l}x_{l,1}}}} + {\sum\limits_{{k \in S_{ny}},{k \neq 0},{l \in S_{nx}}}\;{c_{2{kl}}y_{k,1}x_{l,1}}} + {y_{0}{\sum\limits_{l,{i \in S_{nx}}}\;{c_{30{li}}x_{l,1}x_{i,1}}}} + {\sum\limits_{{k \in S_{ny}},{k \neq 0},l,{i \in S_{nx}}}\;{c_{3{kli}}y_{k,1}x_{l,1}x_{i,1}}}} = {{\sum\limits_{l \in S_{nx}}\;{b_{1l}x_{l,1}^{3}}} + {\sum\limits_{l,{k \in S_{nx}},{l \neq k}}\;{b_{2{lk}}x_{l,1}^{2}x_{k,1}}} + {\sum\limits_{l,{k \in S_{nx}}}\;{b_{3{lk}}x_{l,1}x_{k,1}}} + {\sum\limits_{k \in S_{nx}}\;{b_{4k}x_{k,1}}} + b_{c}}},$

-   -   -   which corresponds to point inverse; and

    -   iii. all other conditions except those in Cases B.i and B.ii.         This case occurs only when ny is greater than or equal to one.

For Case A, for at least one jεS_(nx) x_(j,1)≠x_(j,2), a straight line in (nx+ny+2)-dimensional space is defined by

${\frac{y_{k} - y_{k,1}}{y_{k,2} - y_{k,1}} = \frac{x_{j} - x_{j,1}}{x_{j,2} - x_{j,1}}},{k \in S_{ny}}$ and jεS_(nx) and

${\frac{x_{i} - x_{i,1}}{x_{i,2} - x_{i,1}} = \frac{x_{j} - x_{j,1}}{x_{j,2} - x_{j,1}}},{i \neq j},{i \in {S_{nx}.}}$

For this case, y_(k)=m_(yk)x_(j)+c_(yk), where

$m_{yk} = \frac{y_{k,2} - y_{k,1}}{x_{j,2} - x_{j,1}}$ and c_(yk)=y_(k,1)−x_(j,1)m_(yk). Further, x_(i)=m_(xi)x_(j)+c_(xi), where

$m_{xi} = \frac{x_{i,2} - x_{i,1}}{x_{j,2} - x_{j,1}}$ and c_(xi)=x_(i,1)−x_(j,1)m_(xi). Equation (1) can then be re-written as:

${{{\sum\limits_{k \in S_{ny}}\;{a_{1k}y_{k}^{2}}} + {\sum\limits_{k,{l \in S_{ny}},{l \neq k}}\;{a_{2{kl}}y_{k}y_{l}}} + {\sum\limits_{k \in S_{ny}}\;{a_{3k}y_{k}}} + {x_{j}{\sum\limits_{k,{l \in S_{ny}}}\;{c_{1{klj}}y_{k}y_{l}}}} + {\sum\limits_{k,{l \in S_{ny}},{i \in S_{nx}},{i \neq j}}\;{c_{1{kli}}y_{k}y_{l}x_{i}}} + {x_{j}{\sum\limits_{k \in S_{ny}}\;{c_{2{kj}}y_{k}}}} + {\sum\limits_{{k \in S_{ny}},{l \in S_{nx}},{l \neq j}}\;{c_{2{kl}}y_{k}x_{l}}} + {x_{j}^{2}{\sum\limits_{k \in S_{ny}}\;{c_{3{kjj}}y_{k}}}} + {x_{j}{\sum\limits_{{k \in S_{ny}},{l \in S_{nx}},{l \neq j}}\;{c_{3{klj}}y_{k}x_{l}}}} + {x_{j}{\sum\limits_{{k \in S_{ny}},{i \in S_{nx}},{i \neq j}}\;{c_{3{kji}}y_{k}x_{i}}}} + {\sum\limits_{{k \in S_{ny}},l,{i \in S_{nx}},{{{l\&}i} \neq j}}\;{c_{3{kli}}y_{k}x_{l}x_{i}}}} = {{b_{1j}x_{j}^{3}} + {\sum\limits_{{l \in S_{nx}},{l \neq j}}\;{b_{1l}x_{l}^{3}}} + {x_{j}^{2}{\sum\limits_{{k \in S_{nx}},{k \neq j}}\;{b_{2{jk}}x_{k}}}} + {x_{j}{\sum\limits_{{l \in S_{nx}},{l \neq j}}\;{b_{2{lj}}x_{l}^{2}}}} + {\sum\limits_{l,{k \in S_{nx}},l,{k \neq j},{l \neq k}}\;{b_{2{lk}}x_{l}^{2}x_{k}}} + {b_{3{jj}}x_{j}^{2}} + {x_{j}{\sum\limits_{{k \in S_{nx}},{k \neq j}}\;{b_{3{jk}}x_{k}}}} + {x_{j}{\sum\limits_{{l \in S_{nx}},{l \neq j}}\;{b_{3{lj}}x_{l}}}} + {\sum\limits_{l,{k \in S_{nx}},l,{k \neq j}}\;{b_{3{lk}}x_{l}x_{k}}} + {b_{4j}x_{j}} + {\sum\limits_{{k \in S_{nx}},{k \neq j}}\;{b_{4k}x_{k}}} + b_{c}}},$ and substitution of the above into the rewritten equation (1) for y_(k), kεS_(ny) and x_(i), iεS_(nx) & i≠j, results in:

${{\sum\limits_{k \in S_{ny}}\;{a_{1k}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}^{2}} + {\sum\limits_{k,{l \in S_{ny}},{l \neq k}}\;{{a_{2{kl}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}\left( {{m_{yl}x_{j}} + c_{yl}} \right)}} + {\sum\limits_{k \in S_{ny}}\;{a_{3k}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}} + {x_{j}{\sum\limits_{k,{l \in S_{ny}}}\;{{c_{1{klj}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}\left( {{m_{yl}x_{j}} + c_{yl}} \right)}}} + {\sum\limits_{k,{l \in S_{ny}},{i \in S_{nx}},{i \neq j}}\;{{c_{1{kli}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}\left( {{m_{yl}x_{j}} + c_{yl}} \right)\left( {{m_{xi}x_{j}} + c_{xi}} \right)}} + {x_{j}{\sum\limits_{k \in S_{ny}}\;{c_{2{kj}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}}} + {\sum\limits_{{k \in S_{ny}},{l \in S_{nx}},{l \neq j}}\;{{c_{2{kl}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}\left( {{m_{xl}x_{j}} + c_{xl}} \right)}} + {x_{j}^{2}{\sum\limits_{k \in S_{ny}}\;{c_{3{kjj}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}}} + {x_{j}{\sum\limits_{{k \in S_{ny}},{l \in S_{nx}},{l \neq j}}\;{{c_{3{klj}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}\left( {{m_{xl}x_{j}} + c_{xl}} \right)}}} + {x_{j}{\sum\limits_{{k \in S_{ny}},{i \in S_{nx}},{i \neq j}}\;{{c_{3{kji}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}\left( {{m_{xi}x_{j}} + c_{xi}} \right)}}} + {\sum\limits_{{k \in S_{ny}},l,{i \in S_{nx}},{{{l\&}i} \neq j}}\;{{c_{3{kli}}\left( {{m_{yk}x_{j}} + c_{yk}} \right)}\left( {{m_{xl}x_{j}} + c_{xl}} \right)\left( {{m_{xi}x_{j}} + c_{xi}} \right)}}} = {{b_{1j}x_{j}^{3}} + {\sum\limits_{{l \in S_{nx}},{l \neq j}}\;{b_{1l}\left( {{m_{xl}x_{j}} + c_{l}} \right)}^{3}} + {x_{j}^{2}{\sum\limits_{{k \in S_{nx}},{k \neq j}}\;{b_{2{jk}}\left( {{m_{xk}x_{j}} + c_{xk}} \right)}}} + {x_{j}{\sum\limits_{{l \in S_{nx}},{l \neq j}}\;{b_{2{lj}}\left( {{m_{xl}x_{j}} + c_{xl}} \right)}^{2}}} + {\sum\limits_{l,{k \in S_{nx}},{{{l\&}k} \neq j},{l \neq k}}\;{{b_{2{lk}}\left( {{m_{xl}x_{j}} + c_{xl}} \right)}^{2}\left( {{m_{xk}x_{j}} + c_{xk}} \right)}} + {b_{3{jj}}x_{j}^{2}} + {x_{j}{\sum\limits_{{k \in S_{nx}},{k \neq j}}\;{b_{3{jk}}\left( {{m_{xk}x_{j}} + c_{xk}} \right)}}} + {x_{j}{\sum\limits_{l,{k \in S_{nx}},{l \neq j}}\;{b_{3{lj}}\left( {{m_{xl}x_{j}} + c_{xl}} \right)}}} + {\sum\limits_{l,{k \in S_{nx}},{{{l\&}k} \neq j}}\;{{b_{3{lk}}\left( {{m_{xl}x_{j}} + c_{xl}} \right)}\left( {{m_{xk}x_{j}} + c_{xk}} \right)}} + {b_{4j}x_{j +}{\sum\limits_{{k \in S_{nx}},{k \neq j}}\;{a_{6k}\left( {{m_{xk}x_{j}} + c_{xk}} \right)}}} + b_{c}}$

Expanding the terms in the above equation leads to a cubic equation of the form x_(j), c₃x_(j) ³+C₂x_(j) ²+C₁x_(j)+C₀=0, where C₃, c₂, c₁ & C₀ are obtained from the above equation.

Assuming C₃≠0, the above cubic equation in x_(j) has three roots x_(j,1)x_(j,2), & x_(j,3)′ and can be written as (x_(j)−x_(j,1))(x_(j)−x_(j,2))(x_(j)−x_(j,3)′)=0. Normalizing by the coefficient of x³ and equating the coefficients of x² in the resulting equation with that in (X_(j)−x_(j,1))(x_(j)−x_(j,2))(x_(j)−x_(j,3)′)=0, one obtains a solution for x_(j,3)′:

$\begin{matrix} {x_{j,3}^{\prime} = {\frac{- C_{2}}{C_{3}} - x_{j,1} - {x_{j,2}.}}} & (2) \end{matrix}$ The values of y_(k,3)′,kεS_(ny), and x_(i,3),iεS_(nx) & i≠j, may be similarly obtained from equations for x_(j)=x_(j,3)′.

For cases where C₃=0, c₃x_(j) ³+C₂x_(j) ²+C₁x_(j)+C₀=0 becomes a quadratic equation. Such quadratic equations may be used in the definition of point equivalences.

With regard to Case B for all jεS_(nx),x_(j,1)=x_(j,2), the three sub-cases are considered below. In all cases, x_(j,0) is defined as x_(j,0)=x_(j,1)=x_(j,2), jεS_(nx).

For Case B.i., all kεS_(ny),y_(k,1)=y_(k,2), which corresponds to point doubling. In this case, (x_(0,1), x_(1,1), . . . , x_(nx,1), y_(0,1), y_(1,1), . . . , y_(ny,1))=(x_(0,2), x_(1,2), . . . , x_(nx,2), y_(0,2), y_(1,2), . . . , y_(ny,2)) Letting: (x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,o) ,y _(1,o) , . . . , y _(ny,o))=(x _(0,1) ,x _(1,1) , . . . , x _(nx,1) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))=(x _(0,2) ,x _(1,2) , . . . , x _(nx,2) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2)) the sum is written as (x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,o) ,y _(1,o) , . . . , y _(ny,o))+(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,o) ,y _(1,o) , . . . , y _(ny,o))  (3).

There are several ways of defining the addition in this case. Three possible rules are described below. Case B.i.1: Letting S_(nx,Lx) denote a subset of S_(nx) with Lx elements, i.e., S_(nx,Lx) ⊂S_(nx); letting S_(ny,Ly) denote a subset of S_(ny) with Ly elements and which does not include the element 0; i.e. S_(ny,Ly) ⊂S_(ny) and 1∉S_(ny,Ly); setting the value of Lx and Ly as at least one, then the straight line in this case can be defined as a tangent to the point (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,o), y_(1,o), . . . , y_(ny,o)) defined in a sub-dimensional space with coordinates y_(n) and x_(m) with nεS_(ny,Ly) and mεS_(nx,Lx).

In this case, the gradients m_(yn) and m_(xm) n of the straight line to be used in equation (2) are essentially the first derivatives of y_(n) and x_(m), nεS_(ny,Ly) and mεS_(nx,Lx), for F with respect to x_(j), jεS_(nx,Lx) i.e.,

$m_{yn} = {{\frac{\mathbb{d}y_{n}}{\mathbb{d}x_{j}}\mspace{14mu}{and}\mspace{14mu} m_{xn}} = {\frac{\mathbb{d}x_{m}}{\mathbb{d}x_{j}}.}}$

Using these derivatives for the values of the gradients,

${m_{yn} = \frac{\mathbb{d}y_{n}}{\mathbb{d}x_{j}}},$ where nεS_(ny,Ly), and

${m_{xn} = \frac{\mathbb{d}x_{m}}{\mathbb{d}x_{j}}},$ where MεS_(nx,Lx), in equation (2) and noting that it is assumed that

${\frac{\mathbb{d}x_{m}}{\mathbb{d}x_{j}} = 0},{{{for}\mspace{14mu} m} \in \left( {S_{nx} - S_{{nx},{Lx}}} \right)}$ and ${\frac{\mathbb{d}y_{n}}{\mathbb{d}x_{j}} = 0},{{{for}\mspace{14mu} n} \in \left( {S_{ny} - S_{{ny},{Lx}}} \right)},$ then a solution for x_(j,3)′ may be obtained.

The values of y_(n,3)′ for nεS_(ny) and x_(m,3)′ for mεS_(nx) & m≠j, can further be obtained for x_(j)=x_(j,3)′. The choice of the x_(m)-coordinates, mεS_(nx,Lx), and y_(n)-coordinates, nεS_(ny,Ly), which can be used to compute the tangent of the straight line in this case may be chosen at random or according to a pre-defined rule. Further, a different choice of the x_(m)-coordinates, mεS_(nx,Lx) and y_(n)-coordinates, nεS_(ny,Ly), may be made when one needs to compute successive point doublings, such as that needed in scalar multiplication.

With regard to the next case, Case B.i.2, the second possible way of defining the addition of a point with itself is to apply a sequence of the point doublings according to the rule defined above in Case B.i.1, where the rule is applied with a different selection of the x-coordinate(s) and y-coordinates(s) in each step of this sequence.

In the third sub-case, Case B.i.3, a point is substituted with one of its equivalents. Letting (x_(0,oe), x_(1,oe), . . . , x_(nx,oe), y_(0,oe), y_(1,oe), . . . , y_(ny,oe)) represent the equivalent point of (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,o), y_(1,o), . . . , y_(ny,o)), then equation (3) may be written as: (x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,o) ,y _(1,o) , . . . , y _(ny,o))+(x _(0,oe) ,x _(1,oe) , . . . , x _(nx,oe) ,y _(0,oe) ,y _(1,oe) , . . . , y _(ny,oe)).

With regard to Case B.ii, for kεS_(ny) & k≠0,y_(k,1)=y_(k,2), and where y_(0,1) & y_(0,2) are the roots of the quadratic equation in y₀, this case corresponds to generation of the point inverse.

Letting y_(k,1)=y_(k,2)=y_(k,o) for kεS_(ny) & k≠0, then any two points, such as the point (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,1), y_(1,o), . . . , y_(ny,o))εEC^(nx+ny+2) and the point (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,2), y_(1,o), . . . , y_(ny,o))εEC^(nx+ny+2), are in the hyper-plane with x_(i)=x_(i,o),iεS_(nx) and y_(k)=y_(k,o), kεS_(ny) & k≠0. Thus, any straight line joining these two points such that (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,1), y_(1,o), . . . , y_(ny,o))≠(x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,2), y_(1,o), . . . , y_(ny,o)) is also in this hyper-plane.

Substituting the values of x_(0,o), x_(1,o), . . . , x_(nx,o), y_(1,o), . . . , & y_(ny,o) in an elliptic polynomial equation with multiple x-coordinates and multiple y-coordinates, a quadratic equation for y₀ is obtained, as given above. Thus, y₀ has only two solutions, y_(0,1), & y_(0,2).

Thus, a line joining points (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,1), y_(1,o), . . . , y_(ny,o))εEC^(nx+ny+2) and (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,2), y_(1,o), . . . , y_(ny,o))εEC^(nx+ny+2) does not intersect with EC^(nx+ny+2) at a third point.

A line that joins these two points is assumed to intersect with EC^(nx+ny+2) at infinity (x_(0,I), x_(1,I), . . . , x_(nx,I), y_(0,I), y_(1,I), . . . , y_(ny,I))εEC^(nx+ny+2). This point at infinity is used to define both the inverse of a point in EC^(nx+ny+2) and the identity point. According to the addition rule defined above, one can write: (x ₀ ,x ₁ , . . . , x _(nx) ,y _(0,1) y ₁ , . . . , y _(ny))+(x ₀ ,x ₁ , . . . , x _(nx) ,y _(0,2) ,y ₁ , . . . , y _(ny))=(x _(0,I) ,x _(1,I) , . . . , x _(nx,I) ,y _(0,I) ,y _(1,I) , . . . , y _(ny,I))  (4), since the third point of intersection of such lines is assumed to be at infinity, (x_(0,I), x_(1,I), . . . , x_(nx,I), y_(0,I), y_(1,I), . . . , y_(ny,I))εEC^(nx+ny+2). Thus, this equation defines a unique inverse for any point (x₀, x₁, . . . , x_(nx), y₀, y₁, . . . , y_(ny))εEC^(nx+ny+2), namely: −(x ₀ , x ₁ , . . . , x _(nx) , y _(0,1) , y ₁ , . . . , y _(ny))=(x ₀ , x ₁ , . . . , x _(nx) , y _(0,2) , y ₁ , . . . , y _(ny)).

Thus, equation (4) can be written as: (x ₀ ,x ₁ , . . . , x _(nx) ,y _(0,1) ,y ₁ , . . . , y _(ny))−(x ₀ ,x ₁ , . . . , x _(nx) ,y _(0,1) ,y ₁ , . . . , y _(ny))=(x _(0,I) ,x _(1,I) , . . . , x _(nx,I) ,y _(0,I) ,y _(1,I) , . . . , y _(ny,I))  (5).

Further, a line joining the point at infinity (x_(0,I), x_(1,I), . . . , x_(nx,I), y_(0,I), y_(1,I), . . . , y_(ny,I))εEC^(nx+ny+2) and a point (x₀, x₁, . . . , x_(nx), y_(0,1), y₁, . . . , y_(ny))εEC^(nx+ny+2) will intersect with EC^(nx+ny+2) at (x₀, x₁, . . . , x_(nx), y_(0,2), y₁, . . . , y_(ny))εEC^(nx+ny+2). Thus, from the addition rule defined above, (x ₀ ,x ₁ , . . . , x _(nx) ,y ₀ ,y ₁ ,y ₂ , . . . , y _(ny))+(x _(0,I) ,x _(1,I) , . . . , x _(nx,I) ,y _(0,I) ,y _(1,I) , . . . , y _(ny,I))=(x ₀ ,x ₁ , . . . , x _(nx) ,y ₀ ,y ₁ , . . . , y _(ny))  (6). Equation (5) satisfies axiom (ii) while equation (6) satisfies axiom (i), defined above.

Case B.iii applies for all other conditions except those in cases B.i and B.ii. This case occurs only when ny is greater than or equal to one. Given two points (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,1), y_(1,1), . . . , y_(ny,1))εEC^(nx+ny+2) and (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,2), y_(1,2), . . . , y_(ny,2))εEC^(nx+ny+2) that do not satisfy the conditions of cases B.i and B.ii above, the sum point is written as (x_(0,3), x_(1,3), . . . , x_(nx,3), y_(0,3), y_(1,3), . . . , y_(ny,3))=(x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,1), y_(1,1), . . . , y_(ny,1))+(x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,2), y_(1,2), . . . , y_(ny,2)).

There are several possible rules to find the sum point in this case. Three possible methods are given below.

1) Using three point doublings and one point addition, (x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=4(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))−2(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2));

2) using one point doublings and three point additions,

(x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=(2(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))+(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,2) ,y _(1,2) , . . . , y _(ny,2)))−(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1)); and

3) using point equivalence, (x _(0,3) ,x _(1,3) , . . . , x _(nx,3) ,y _(0,3) ,y _(1,3) , . . . , y _(ny,3))=(x _(0,o) ,x _(1,o) , . . . , x _(nx,o) ,y _(0,1) ,y _(1,1) , . . . , y _(ny,1))+(x _(0,oe) ,x _(1,oe) , . . . , x _(nx,oe) ,y _(0,2e) ,y _(1,2e) , . . . , y _(ny,2e)), where (x_(0,oe), x_(1,oe), . . . , x_(nx,oe), y_(0,2e), y_(1,2e), . . . , y_(ny,2e)) is assumed to be the equivalent point of (x_(0,o), x_(1,o), . . . , x_(nx,o), y_(0,2), y_(1,2), . . . , y_(ny,2)).

It should be noted that the above methods for defining the sum point are not the only ones that can be defined and are provided for exemplary purposes only. The choice of method used to obtain the sum point in this case should depend on the computation complexity of point addition and point doubling.

With regard to associativity, one way of proving associativity of (EC^(nx+,ny+2),+) is as follows: Given particular elliptic polynomial equations (i.e., for particular coefficients a_(1l), a_(2kl), a_(3k), c_(1lki), c_(2kl), c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k), b_(c)εF) defined over a finite field F, if it can be shown by computation that any point QεEC^(nx+ny+2) (and any of its equivalent points) can be uniquely written as k_(Q)PεEC^(nx+ny+2), where P is the generator point of (EC^(nx+ny+2),+), then the corresponding EC^(nx+ny+2) groups based on such polynomials are associative. This is because any-three points Q, R, SεEC^(nx+ny+2) (or any of their equivalent points) can be written as k_(Q)P, k_(R)P, k_(S)PεEC^(nx+ny+2), respectively, thus their sum (Q+R+S)=(k_(Q)P+k_(R)P+k_(S)P)=(k_(Q)+k_(R)+k_(S)), can be carried out in any order.

The following elliptic polynomial equation with nx=1 and ny=0 is used to show an example of the equations in Case A used in point addition: y₀ ²=x₀ ³+x₁ ³+x₀x₁. Choosing x_(j)=x₀, and substituting y_(k)=m_(yk)x_(j)+c_(yk) from Case A above for y₀, and the corresponding equation x_(i)=m_(xi)x_(j)+c_(xi) for x₁, one obtains (m_(y0)x₀+c_(y0))²=x₀ ³+(m_(x1)x₀+c_(x1))³+x₀(m_(x1)x₀+c_(x1)).

Expanding this equation yields the equation m_(y0) ²x₀ ²+2m_(y0)c_(y0)x₀+c_(y0) ²=x₀ ³+m_(x1) ³x₀ ³+3m_(x1) ²c_(x1) ²+3m_(x1)c_(x1) ²x₀+c_(x1) ³+m_(x1)x₀ ²+c_(x1)x₀, or (1+m_(x1) ³)x₀ ³+(3m_(x1) ^(2c) _(x1)+m_(x1)−m_(y0) ²)x₀ ²+(3m_(x1)c_(x1) ²+c_(x1)−2m_(y0)c_(y0))x₀+c_(x1) ³−c_(y0) ²=0. From equation (2), the solution for x_(0,3)′ in this case is obtained:

$x_{0,3}^{\prime} = {\frac{- \left( {{3m_{x\; 1}^{2}c_{x\; 1}} + m_{x\; 1} - m_{y\; 0}^{2}} \right)}{\left( {1 + m_{x\; 1}^{3}} \right)} - x_{j,1} - {x_{j,2}.}}$ Similarly, one can obtain the values of y_(0,3) and x_(1,3) for X₀=x_(0,3).

It should be noted that when m_(x1)=−1, the coefficient of the cubic term in the above is zero; i.e. C₃=0. In this case, the resulting quadratic equation can be used in the definition of point equivalences for the points that satisfy the elliptic polynomial equation.

Each of the equations for point addition and point doublings derived for cases A and B above require modular inversion or division. In cases where field inversions or divisions are significantly more expensive (in terms of computational time and energy) than multiplication, projective coordinates are used to remove the requirement for field inversion or division from these equations.

Several projective coordinates can be utilized. In the preferred embodiment, the Jacobean projective coordinate system is used. As an example, we examine:

$\begin{matrix} {{x_{i} = \frac{X_{i}}{V^{2}}}{{{{for}\mspace{14mu} i} \in S_{nx}};}{{and};}{and}} & (7) \\ {y_{k} = {{\frac{Y_{k}}{V^{3}}\mspace{14mu}{for}\mspace{14mu} k} \in {S_{ny}.}}} & (8) \end{matrix}$

Using Jacobian projection yields:

$\begin{matrix} {{{\sum\limits_{k \in S_{ny}}\;{a_{1k}\frac{Y_{k}^{2}}{V^{6}}}} + {\sum\limits_{k,{l \in S_{ny}},{l \neq k}}\;{a_{2k\; l}\frac{Y_{k}}{V^{3}}\frac{Y_{l}}{V^{3}}}} + {\sum\limits_{k \in S_{ny}}\;{a_{3k}\frac{Y_{k}}{V^{3}}}} + {\sum\limits_{k,{l \in S_{ny}},{i \in S_{nx}}}\;{c_{1{kli}}\frac{Y_{k}}{V^{3}}\frac{Y_{l}}{V^{3}}\frac{X_{i}}{V^{2}}}} + {\sum\limits_{{k \in S_{ny}},{l \in S_{nx}}}\;{c_{2{kl}}\frac{Y_{k}}{V^{3}}\frac{X_{l}}{V^{2}}}} + {\sum\limits_{{k \in S_{ny}},l,{i \in S_{nx}}}\;{c_{3{kli}}\frac{Y_{k}}{V^{3}}\frac{X_{l}}{V^{2}}\frac{X_{i}}{V^{2}}}}} = {{\sum\limits_{l \in S_{nx}}\;{b_{1l}\frac{X_{l}^{3}}{V^{6}}}} + {\sum\limits_{l,{k \in S_{nx}},{l \neq k}}\;{b_{2{lk}}\frac{X_{l}^{2}}{V^{4}}\frac{X_{k}}{V^{2}}}} + {\sum\limits_{l,{k \in S_{nx}}}\;{b_{3{lk}}\frac{X_{l}}{V^{2}}\frac{X_{k}}{V^{2}}}} + {\sum\limits_{k \in S_{nx}}\;{b_{4k}\frac{X_{k}}{V^{2}}}} + b_{c}}} & (9) \end{matrix}$ which can be rewritten as:

$\begin{matrix} {{{\sum\limits_{k \in S_{ny}}\;{a_{1k}Y_{k}^{2}V^{2}}} + {\sum\limits_{k,{l \in S_{ny}},{l \neq k}}\;{a_{2{kl}}Y_{k}Y_{l}V^{2}}} + {\sum\limits_{k \in S_{ny}}\;{a_{3k}Y_{k}V^{5}}} + {\sum\limits_{k,{l \in S_{ny}},{i \in S_{nx}}}\;{c_{1{kli}}Y_{k}Y_{l}X_{i}}} + {\sum\limits_{{k \in S_{ny}},{l \in S_{nx}}}\;{c_{2{kl}}Y_{k}X_{l}V^{3}}} + {\sum\limits_{{k \in S_{ny}},l,{i \in S_{nx}}}\;{c_{3{kli}}Y_{k}X_{l}X_{i}V}}} = {{\sum\limits_{l \in S_{nx}}\;{b_{1l}X_{l}^{3}V^{2}}} + {\sum\limits_{l,{k \in S_{nx}},{l \neq k}}\;{b_{2{lk}}X_{l}^{2}X_{k}V^{2}}} + {\sum\limits_{l,{k \in S_{nx}}}\;{b_{3{lk}}X_{l}X_{k}V^{4}}} + {\sum\limits_{k \in S_{nx}}\;{b_{4k}X_{k}V^{6}}} + {b_{c}{V^{8}.}}}} & (10) \end{matrix}$

In the following, the points (X₀, X₁, . . . , X_(nx), Y₀, Y₁, . . . , Y_(ny), V) are assumed to satisfy equation (10). When V≠0, the projected point (X₀, X₁, . . . , X_(nx), Y₀, Y₁, . . . , Y_(ny), V) corresponds to the point:

${\left( {x_{0},x_{1},\ldots\mspace{11mu},x_{nx},y_{0},y_{1},\ldots\mspace{11mu},y_{ny}} \right) = \left( {\frac{X_{0}}{V^{2}},\frac{X_{1}}{V^{2}},\ldots\mspace{11mu},\frac{X_{nx}}{V^{2}},\frac{Y_{0}}{V^{3}},\frac{Y_{1}}{V^{3}},\ldots\mspace{11mu},\frac{Y_{ny}}{V^{3}}} \right)},$ which satisfies equation (1).

Using Jacobean projective coordinates, equation (10) can be written as:

$\begin{matrix} {\left( {\frac{X_{0,3}}{V_{3}^{2}},\frac{X_{1,3}}{V_{3}^{2}},\ldots\mspace{11mu},\frac{X_{{nx},3}}{V_{3}^{2}},\frac{Y_{0,3}}{V_{3}^{3}},\frac{Y_{1,3}}{V_{3}^{3}},\ldots\mspace{11mu},\frac{Y_{{ny},3}}{V_{3}^{3}}} \right) = {\left( {\frac{X_{0,1}}{V_{1}^{2}},\frac{X_{1,1}}{V_{1}^{2}},\ldots\mspace{11mu},\frac{X_{{nx},1}}{V_{1}^{2}},\frac{Y_{0,1}}{V_{1}^{3}},\frac{Y_{1,1}}{V_{1}^{3}},\ldots\mspace{11mu},\frac{Y_{{ny},1}}{V_{1}^{3}}} \right) + {\left( {\frac{X_{0,2}}{V_{2}^{2}},\frac{X_{1,2}}{V_{2}^{2}},\ldots\mspace{11mu},\frac{X_{{nx},2}}{V_{2}^{2}},\frac{Y_{0,2}}{V_{2}^{3}},\frac{Y_{1,2}}{V_{2}^{3}},\ldots\mspace{11mu},\frac{Y_{{ny},2}}{V_{2}^{3}}} \right).}}} & (11) \end{matrix}$

By using Jacobian projective coordinates in the equations of Cases A and B above, and by an appropriate choice of the value of V₃, it can be shown that point doubling and point addition can be computed without the need for field inversion or division.

One example of an elliptic polynomial and its twist is given below. Writing the cubic equation with two variables in x-coordinates and one variable in y-coordinates defined over a finite field F yields: t=x ₁ ³ +ax ₂ ³ +bx ₁ x ₂ +c  (12), where x₁ & x₂εF tεF and a, b & cεF. Any value of x₁,x₂ will lead to a value of tεF. It should be noted that t could be quadratic residue or non-quadratic residue. If t is quadratic residue, it can be written as t=y², and if t is non-quadratic residue, it can be written as t= αy², where α is a non-quadratic element of F; i.e., √{square root over ( α)}∉F. Therefore, equation (12) can be written as:

$\begin{matrix} {{{\alpha\; y^{2}} = {x_{1}^{3} + {ax}_{2}^{3} + {{bx}_{1}x_{2}} + c}}{where}\left\{ {\begin{matrix} {\alpha = 1} & {{if}\mspace{14mu} t\mspace{14mu}{is}\mspace{14mu}{quadraticresidue}} \\ {\alpha = \overset{\_}{\alpha}} & {{if}\mspace{14mu} t\mspace{14mu}{is}\mspace{14mu}{non}\text{-}{quadraticresidue}} \end{matrix}.} \right.} & (13) \end{matrix}$

It should be noted that for specific coefficient a, b & cεF, when α=1, the resulting polynomial is an elliptic polynomial. However, if α= α this leads to a twist of the elliptic polynomial obtained with α=1.

Any value of x₁ & x₂ εF will lead to a point (x₁,x₂,√{square root over (α)}y), which is either on an elliptic polynomial or its twist. If α=1, the point (x₁,x₂,√{square root over (α)}y) is on the elliptic polynomial. If α== α, the point (x₁,x₂,√ αy) is on its twist.

Elliptic point addition can be formulated on a twist of an elliptic polynomial in the same fashion as it is formulated for elliptic polynomials and described above. As a result, elliptic polynomial cryptography can also be defined on twists of elliptic polynomial in the same fashion as that described above.

In the following MAC generation methods, a shared secret key or keys are used to generate the MAC of the message data bits.

In a first embodiment, the MAC generation method includes the following steps:

a) defining a maximum block size that can be embedded into (nx+1) x-coordinates and ny y-coordinates, wherein nx and ny are integers, and setting the maximum block size to be (nx+ny+1)N bits, wherein N is an integer, and defining the number of (nx+ny+1)N-bit blocks in a message bit string to be (u+1), wherein u is an integer;

b) a sending correspondent and a receiving correspondent agree upon the values of nx and ny, and further agree on a set of coefficients a_(1k), a_(2kl), a_(3k), c_(1lki), c_(2kl), c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k) & b_(c)εF, wherein F represents a finite field where the field's elements can be represented in N-bits, the sending and receiving correspondents further agreeing on a random number k that represents a shared secret key for communication, the sending and receiving correspondents further agreeing on a base point on an elliptic polynomial (x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B))εEC^(nx+ny+2) wherein EC^(nx+ny+2) represents the elliptic polynomial, and further agreeing on a base point on the twist of the elliptic polynomial (x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB))εTEC^(xn+ny+2), wherein TEC^(nx+ny+2) represents the twist of the elliptic polynomial EC^(nx+ny+2), and further, α is a non-quadratic residue element of F;

the sending correspondent then performs the following steps:

c) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of a key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using a non-iterative embedding method to be described below;

d) if (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, a scalar multiplication is computed as the equations: (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=k(x _(0,TB) ,x _(1,TB) , . . . , x _(nx,TB) ,√{square root over ( α)}y _(0,B) ,y _(1,TB) , . . . , y _(ny,TB)) and (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,y _(0,k) ,y _(1,k) , . . . , y _(ny,k)). Otherwise, if α_(k)= α, then the scalar multiplication is computed as: (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=k(x _(0,B) ,x _(1,B) , . . . , x _(nx,B) ,y _(0,B) ,y _(1,B) , . . . , y _(ny,B)) and (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,√{square root over ( α)}y _(0,k) ,y _(1,k) , . . . , y _(ny,k));

e) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), where M and L are integers, and embedding the message bit string of the initial block of the message bit string into the (nx+1) x-coordinates and the ny y-coordinates of the elliptic message point (x_(0,m) ₀ , x_(1,m) ₀ , . . . , x_(nx,m) ₀ , √{square root over (α)}y_(0,m) ₀ , y_(1,m) ₀ , . . . , y_(ny,m) ₀ ) using the non-iterative embedding method;

f) if the message point of the 0-th block is on the elliptic polynomial, then the message authentication code points are computed as (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,√{square root over ( α)}y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,√{square root over ( α)}y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,√{square root over ( α)}y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x_(0,Tc) ₀ , x_(1,Tc) ₀ , . . . , x_(nx,Tc) ₀ , √{square root over ( α)}y_(0,Tc) ₀ , y_(1,Tc) ₀ , . . . , y_(ny,Tc) ₀ )=(x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ ). Otherwise, they are computed as the equations: (x _(0,Tc) ₀ ,x _(1,Tc) ₀ , . . . , x _(nx,Tc) ₀ ,y _(0,Tc) ₀ ,y _(1,Tc) ₀ , . . . , y _(ny,Tc) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ ) and (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ );

g) initializing an integer i as i=1, and repeating the following steps h) and i) until i>u:

h) embedding the message (nx+ny+1)N-bit string of an i-th block into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx) and the ny y-coordinates, y₁, . . . , y_(ny), of the elliptic message point (x_(0,m) _(i) , x_(1,m) _(i) , . . . , x_(nx,m) _(i) , √{square root over (α_(m) _(i) )}y_(0,m) _(i) , y_(1,m) _(i) , . . . , y_(ny,m) _(i) ) using the non-iterative embedding method

i) if the message point of the i-th block is on the elliptic polynomial, then the message authentication code points are computed as: (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i) , . . . , x _(ny,c) _(i−1) ) and (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,√{square root over ( α)}y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , x _(ny,Tc) _(i−1) ). Otherwise, they are computed as (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,√{square root over ( α)}y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,√{square root over ( α)}y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , x _(ny,Tc) _(i−1) ) and (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , y _(ny,c) _(i−1) );

j) a set of appropriate bits of the x-coordinates x_(0,c), x_(1,c), . . . , x_(nx,c), y-coordinates y_(1,c), . . . , y_(ny,c) and a sign bit of y_(0,c) of the message authentication code point (x_(0,c) _(u) , x_(1,c) _(u) , . . . , x_(nx,c) _(u) , y_(0,c) _(u) , y_(1,c) _(u) , . . . , y_(ny,c) _(u) ) and a set of appropriate bits of the x-coordinates x_(0,Tc), x_(1,Tc), . . . & x_(nx,Tc), y-coordinates y_(1,Tc), . . . , y_(ny,Tc) and the sign bit of y_(0,Tc) of the message authentication code point (x_(0,Tc) _(u) , x_(1,Tc) _(u) , . . . , x_(nx,Tc) _(u) , y_(0,Tc) _(u) , y_(1,Tc) _(u) , . . . , y_(ny,Tc) _(u) ) are concatenated together to form the message authentication code, which is then appended to the message to be sent together with, if necessary, any additional information needed to help authenticate the message bit string at the receiving correspondent without compromising security and sent to the receiving correspondent;

the receiving correspondent then performs the following steps:

k) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of the key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using the non-iterative embedding method;

l) if the point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, then a scalar multiplication is performed as (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=k(x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB)) and (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), y_(0,k), y_(1,k), . . . , y_(ny,k)) otherwise if α_(k)=α_(o) then the scalar multiplication is performed as (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=k(x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B)) and (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over ( α)}y_(0,k), y_(1,k), . . . , y_(ny,k));

m) dividing the received message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), where M and L are integers, and embedding the received message (nx+ny+1)N-bit string of the initial block of the received message bit string into the (nx+1) x-coordinates and the ny y-coordinates of the received message elliptic point (x_(0,rm) ₀ , x_(1,rm) ₀ , . . . , x_(nx,rm) ₀ , √{square root over (α_(rm) ₀ )}y_(0,rm) ₀ , y_(1,rm) ₀ , . . . , y_(ny,rm) ₀ ) using the non-iterative embedding method;

n) if the message point of the 0-th block is on the elliptic polynomial, then the message authentication points are computed as (x_(0,rc) ₀ , x_(1,rc) ₀ , . . . , x_(nx,rc) ₀ , y_(0,rc) ₀ , y_(1,rc) ₀ , . . . , y_(ny,rc) ₀ )=(x_(0,rm) ₀ , x_(1,rm) ₀ , . . . , x_(nx,rm) ₀ , y_(0,rm) ₀ , y_(1,rm) ₀ , . . . , y_(ny,rm) ₀ )+(x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , y_(ny,S) ₀ ) and (x_(0,rTc) ₀ , x_(1,rTc) ₀ , . . . , x_(nx,rTc) ₀ , √{square root over ( α)}y_(0,rTc) ₀ , y_(1,rTc) ₀ , . . . , y_(ny,rTc) ₀ )=(x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ ). Otherwise, they are computed as: (x _(0,rTc) ₀ ,x _(1,rTc) ₀ , . . . , x _(nx,rTc) ₀ ,√{square root over ( α)}y _(0,rTc) ₀ ,y _(1,rTc) ₀ , . . . , y _(ny,rTc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,√{square root over ( α)}y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,rTS) ₀ ,x _(1,rTS) ₀ , . . . , x _(nx,rTS) ₀ ,√{square root over ( α)}y _(0,rTS) ₀ ,y _(1,rTS) ₀ , . . . , x _(ny,rTS) ₀ ) and (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )

o) initializing the integer i as i=0, the following steps p) and q) are repeated until i>u:

p) embedding the received message (nx+ny+1)N-bit string of the i-th block the into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny) of the received message elliptic point (x_(0,rm) _(i) , x_(1,rm) _(i) , . . . , x_(nx,rm) _(i) , √{square root over (α_(rm) _(i) )}y_(0,rm) _(i) , y_(1,rm) _(i) , . . . , y_(ny,rm) _(i) ) using the non-iterative embedding method;

q) if the received message point of the i-th block is on the elliptic polynomial, then the message authentication points are computed as (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i−1) , . . . , x _(ny,rc) _(i−1) ) and (x_(0,rTc) _(i) , x_(1,rTc) _(i) , . . . , x_(nx,rTc) _(i) , √{square root over ( α)}y_(0,rTc) _(i) , y_(1,rTc) _(i) , . . . , y_(ny,rTc) _(i) )=(x_(0,rTc) _(i−1) , x_(1,rTc) _(i−1) , . . . , x_(nx,rTc) _(i−1) , √{square root over ( α)}y_(0,rTc) _(i−1) , y_(1,rTc) _(i−1) , . . . , x_(ny,rTc) _(i−1) ), otherwise they are computed as (x _(0,rTc) _(i) ,x _(1,rTc) _(i) , . . . , x _(nx,rTc) _(i) ,√{square root over ( α)}y _(0,rTc) _(i) ,y _(1,rTc) _(i) , . . . , y _(ny,rTc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,√{square root over ( α)}y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rTc) _(i−4) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) ,√{square root over ( α)}y _(0,rTc) _(i−1) ,y _(1,rTc) _(i−1) , . . . , x _(ny,rTc) _(i−1) ) and (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i−1) , . . . , y _(ny,rc) _(i−1) )

r) the appropriate bits of the x-coordinates x_(0,rc), x_(1,rc), . . . , x_(nx,rc), y-coordinates y_(1,rc), . . . , y_(ny,rc) and the sign bit of y_(0,rc) of the message authentication code point (x_(0,rc) _(u) , x_(1,rc) _(u) , . . . , x_(nx,rc) _(u) , y_(0,rc) _(u) , y_(1,rc) _(u) , . . . , y_(ny,rc) _(u) ), and the appropriate bits of the x-coordinates x_(0,rTc), x_(1,rTc), . . . , x_(nx,rTC), y-coordinates y_(1,rTc), . . . , y_(ny,rTc) and the sign bit of y_(0,rTc) of the message authentication point (x_(0,rTc) _(u) , x_(1,rTc) _(u) , . . . , x_(nx,rTc) _(u) , √{square root over ( α)}y_(0,rTc) _(u) , y_(1,rTc) _(u) , . . . , y_(ny,rTc) _(u) ) are concatenated together to form the message authentication code of the received message; and

s) if the received message authentication code is equal to the message authentication code of the received message generated at the receiving correspondent, the received message is authenticated.

In an alternative embodiment, the MAC generation method includes the following steps:

a) defining a maximum block size that can be embedded into (nx+1) x-coordinates and ny y-coordinates, wherein n is an integer, and setting the maximum block size to be (nx+ny+1)N bits, wherein N is an integer, and defining the number of (nx+ny+1)N-bit blocks in a message bit string to be (u+1), wherein u is an integer;

b) a sending correspondent and a receiving correspondent agree upon the values of nx and ny, where nx and ny are integers, and further agree on a set of coefficients a_(1k), a_(2kl), a_(3k),c_(1lki),c_(2kl),c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k) & b_(c)εF, wherein F represents a finite field and where the elements of F are represented in N-bits, the sending and receiving correspondents further agreeing on a random number k that represents a shared secret key for communication, the sending and receiving correspondents further agreeing on a base point on an elliptic polynomial (x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B))εEC^(nx+ny+2), wherein EC^(nx+ny+2) represents the elliptic polynomial, and further agreeing on a base point on the twist of the elliptic polynomial (x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB))εTEC^(xn+ny+2), wherein TEC^(nx+ny+2) represents the twist of the elliptic polynomial EC^(nx+ny+2), and where α is a non-quadratic-residue element of F;

the sending correspondent then performs the following steps:

c) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny) of a key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using a non-iterative embedding method to be described below;

d) if (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, a scalar multiplication is computed as: (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=k(x _(0,TB) ,x _(1,TB) , . . . , x _(nx,TB) ,√{square root over ( α)}y _(0,B) ,y _(1,TB) , . . . , y _(ny,TB)) and (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,y _(0,k),y_(1,k) , . . . , y _(ny,k)), otherwise, if α_(k)= α then the scalar multiplication is computed as: (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=(x _(0,B),x_(1,B), . . . , x_(nx,B),y_(0,B),y_(1,B), . . . , y_(ny,B)) and (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )+(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,√{square root over ( α)}y _(0,k) ,y _(1,k) , . . . , y _(ny,k));

e) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), where M and L are integers, and embedding the message bit string of the initial block of the message bit string into the (nx+1) x-coordinates and the ny y-coordinates of the elliptic message point (x_(0,m) ₀ , x_(1,m) ₀ , . . . , x_(nx,m) ₀ , √{square root over (α_(m) ₀ )}y_(0,m) ₀ , y_(1,m) ₀ , . . . , y_(ny,m) ₀ ) using the non-iterative embedding method;

f) if the message point of the 0-th block is on the elliptic polynomial, then the message authentication code points are computed as (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x_(0,Tc) ₀ , x_(1,Tc) ₀ , . . . , x_(nx,Tc) ₀ , √{square root over ( α)}y_(0,Tc) ₀ , y_(1,Tc) ₀ , . . . , y_(ny,Tc) ₀ )=(x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ ). Otherwise, they are computed as: (x _(0,Tc) ₀ ,x _(1,Tc) ₀ , . . . , x _(nx,Tc) ₀ ,√{square root over ( α)}y _(0,Tc) ₀ ,y _(1,Tc) ₀ , . . . , y _(ny,Tc) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,√{square root over ( α)}y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ ) and (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ):

g) initializing an integer i as i=1, and repeating the following steps h) through j) until i>u:

h) embedding the message (nx+ny+1)N-bit string of an i-th block into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates, y₁, . . . , y_(ny), of the elliptic message point (x_(0,m) _(i) , x_(1,m) _(i) , . . . , x_(nx,m) _(i) √{square root over (α_(m) _(i) )}y_(0,m) _(i) , y_(1,m) _(i) , . . . , y_(ny,m) _(i) ) using the non-iterative embedding method;

i) doubling the points (x_(0,c) _(i−1) , x_(1,c) _(i−1) , . . . , x_(nx,c) _(i−1) , y_(0,c) _(i−1) , y_(1,c) _(i−1) , . . . , y_(ny,c) _(i−1) ) and (x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , y_(ny,Tc) _(i−1) ) as: (x_(0,c) _(i−1) , x_(1,c) _(i−1) , . . . , x_(nx,c) _(i−1) , y_(0,c) _(i−1) , y_(1,c) _(i−1) , . . . , y_(ny,c) _(i−1) )=2(x_(0,c) _(i−1) , x_(1,c) _(i−1) , . . . , x_(nx,c) _(i−1) , y_(0,c) _(i−1) , y_(1,c) _(i−1) , . . . , y_(ny,c) _(i−1) ) and (x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , y_(ny,Tc) _(i−1) )=2(x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , y_(ny,Tc) _(i−1) );

j) if the message point of the i-th block is on the elliptic polynomial, then the message authentication code points are computed as (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x_(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , x _(ny,c) _(i−1) ) and (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,√{square root over ( α)}y _(0,Tc) _(i−1) ,y_(1,Tc) _(i−1) , . . . , x_(ny,Tc) _(i−1) ) Otherwise, they are computed as: (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,√{square root over ( α)}y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,√{square root over ( α)}y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , x _(ny,Tc) _(i−1) ) and (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , x _(ny,c) _(i−1) );

k) a set of appropriate bits of the x-coordinates x_(0,c), x_(1,c), . . . , x_(nx,c), y-coordinates y_(1,c), . . . , y_(yn,c) and a sign bit of y_(0,c) of the message authentication code point (x_(0,c) _(u) , x_(1,c) _(u) , . . . , x_(nx,c) _(u) , y_(0,c) _(u) , y_(1,c) _(u) , . . . , y_(ny,c) _(u) ), and a set of appropriate bits of the x-coordinates x_(0,Tc), x_(1,Tc), . . . & x_(nx,Tc) y-coordinates y_(1,Tc), . . . , y_(ny,Tc) and the sign bit of y_(0,Tc) of the message authentication code point (x_(0,Tc) _(u) , x_(1,Tc) _(u) , . . . , x_(nx,Tc) _(u) , √{square root over ( α)}y_(0,Tc) _(u) , y_(1,Tc) _(u) , . . . , y_(ny,Tc) _(u) ) are concatenated together to form the message authentication code, which is then appended to the message to be sent together with, if necessary, any additional information needed to help authenticate the message bit string at the receiving correspondent without compromising security and sent to the receiving correspondent;

the receiving correspondent then performs the following steps:

l) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of the key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using the non-iterative embedding method;

m) if the point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, then a scalar multiplication is performed as (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) _(i) , . . . , x_(ny,TS) ₀ )=k(x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB)) and (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , y_(ny,S) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) otherwise if α_(k)=α_(o) then the scalar multiplication is performed as (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,y _(0,k) ,y _(1,k) , . . . , y _(ny,k)) and (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) _(i) , . . . , x_(ny,TS) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,√{square root over ( α)}y _(0,k) ,y _(1,k) , . . . , y _(ny,k));

n) dividing the received message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), where M and L are integers, and embedding the received message bit string of the initial block of the received message bit string into the (nx+1) x-coordinates and the ny y-coordinates of the received message elliptic point (x_(0,rm) ₀ , x_(1,rm) ₀ , . . . , x_(nx,rm) ₀ , √{square root over (α)}y_(0,rm) ₀ , y_(1,rm) _(i) , . . . , y_(ny,rm) ₀ ) using the non-iterative embedding method;

o) if the message point of the 0-th block is on the elliptic polynomial, then the message authentication points are computed as: (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x_(0,rTc) ₀ , x_(1,rTc) ₀ , . . . , x_(nx,rTc) ₀ , √{square root over ( α)}y_(0,rTc) ₀ , y_(1,rTc) ₀ , . . . , y_(ny,rTc) ₀ )=(x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ ). Otherwise, they are computed as: (x _(0,rTc) ₀ ,x _(1,rTc) ₀ , . . . , x _(nx,rTc) ₀ ,√{square root over ( α)}y _(0,rTc) ₀ ,y _(1,rTc) ₀ , . . . , y _(ny,rTc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,√{square root over ( α)}y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ ) (x _(0,rTS) ₀ ,x _(1,rTS) ₀ , . . . , x _(nx,rTS) ₀ ,√{square root over ( α)}y _(0,rTS) ₀ ,y _(1,rTS) ₀ , . . . , x _(ny,rTS) ₀ ) and (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ , y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )

p) initializing the integer i as i=0, the following steps q) through s) are repeated until i>u:

q) embedding the received message (nx+ny+1)N-bit string of the i-th block the into the (nx+1) x-coordinates x₀,x₁, . . . , x_(nx), and the ny y coordinates y₁, . . . , y_(ny) of the received message elliptic point (x_(0,rm) _(i) , x_(1,rm) _(i) , . . . , x_(nx,rm) _(i) , √{square root over ( α)}y_(0,rm) _(i) , y_(1,rm) _(i) , . . . , y_(ny,rm) _(i) ) using the non-iterative embedding method;

r) doubling the points (x_(0,S) _(i) , x_(1,S) _(i) , . . . , x_(nx,S) _(i) , y_(0,S) _(i) , y_(1,S) _(i) , . . . , x_(ny,S) _(i) ) and (x_(0,TS) _(i) , x_(1,TS) _(i) , . . . , x_(nx,TS) _(i) , √{square root over ( α)}y_(0,TS) _(i) , y_(1,TS) _(i) , . . . , x_(ny,TS) _(i) ) (x _(0,S) _(i) ,x _(1,S) _(i) , . . . , x _(nx,S) _(i) ,y _(0,S) _(i) ,y _(1,S) _(i) , . . . , x _(ny,S) _(i) )=2(x _(0,S) _(i−1) ,x _(1,S) _(i−1) , . . . , x _(nx,S) _(i−1) ,y _(0,S) _(i−1) ,y _(1,S) _(i−1) , . . . , x _(ny,S) _(i−1) ) (x _(0,TS) _(i) ,x _(1,TS) _(i) , . . . , x _(nx,TS) _(i) ,√{square root over ( α)}y _(0,TS) _(i) ,y _(1,TS) _(i) , . . . , x _(ny,TS) ₀ )=2(x _(0,TS) _(i−1) ,x _(1,TS) _(i−1) , . . . , x _(nx,TS) _(i−1) ,√{square root over ( α)}y _(0,TS) _(i−1) ,y _(1,TS) _(i−1) , . . . , x _(ny,TS) _(i−1) );

s) if the received message point of the i-th block is on the elliptic polynomial, then the message authentication points are computed as (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i−1) , . . . , x _(ny,rc) _(i−1) ) and (x_(0,rTc) _(i) , x_(1,rTc) _(i) , . . . , x_(nx,rTc) _(i) , √{square root over ( α)}y_(0,rTc) _(i) , y_(1,rTc) _(i) , . . . , y_(ny,rTc) _(i) )=(x_(0,rTc) _(i−1) , x_(1,rTc) _(i−1) , . . . , x_(nx,rTc) _(i−1) , √{square root over ( α)}y_(0,rTc) _(i−1) , y_(1,rTc) _(i−1) , . . . , x_(ny,rTc) _(i−1) ), otherwise they are computed as: (x _(0,rTc) _(i) ,x _(1,rTc) _(i) , . . . , x _(nx,rTc) _(i) ,√{square root over ( α)}y _(0,rTc) _(i) ,y _(1,rTc) _(i) , . . . , y _(ny,rTc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,√{square root over ( α)}y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rTc) _(i−1) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) ,√{square root over ( α)}y _(0,rTc) _(i−1) ,y _(1,rTc) _(i−1) , . . . , x _(ny,rTc) _(i−1) ) and (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i−1) , . . . , x _(ny,rc) _(i−1) );

t) the appropriate bits of the x-coordinates x_(0,rc), x_(1,rc), . . . , x_(nx,rc), y-coordinates y_(1,rc), . . . , y_(ny,rc) and the sign bit of y_(0,rc) of the message authentication code point (x_(0,rc) _(u) , x_(1,rc) _(u) , . . . , x_(nx,rc) _(u) , y_(0,rc) _(u) , y_(1,rc) _(u) , . . . , y_(ny,rc) _(u) ), and the appropriate bits of the x-coordinates x_(0,rTc), x_(1,rTc), . . . , x_(nx,rTc), Y-coordinates y_(1,rTc), . . . , y_(ny,rTc) and the sign bit of y_(0,rTc) of the message authentication point (x_(0,rTc) _(u) , x_(1,rTc) _(u) , . . . , x_(nx,rTc) _(u) , √{square root over ( α)}y_(0,rTc) _(u) , y_(1,rTc) _(u) , . . . , y_(ny,rTc) _(u) ) are concatenated together to form the message authentication code of the received message; and

u) if the received message authentication code is equal to the message authentication code of the received message generated at the receiving correspondent, the received message is authenticated.

In a further alternative embodiment, the MAC generation method includes the following steps:

a) defining a maximum block size that can be embedded into (nx+1) x-coordinates and ny y-coordinates, wherein n is an integer and nx and ny are integers, and setting the maximum block size to be (nx+ny+1)N bits, wherein N is an integer, and defining the number of (nx+ny+1)N-bit blocks in a message bit string to be (u+1), wherein u is an integer, and further representing the (nx+ny+1)N bit string of the i-th message block, wherein i is an integer, as S_(m) _(i) ;

b) a sending correspondent and a receiving correspondent agree upon the values of nx and ny and further agree on a set of coefficients a_(1k), a_(2kl), a_(3k), c_(1lki), c_(2kl), c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k) & b_(c)εF, wherein F represents a finite field and where the elements of F are represented in N-bits, the sending and receiving correspondents further agreeing on a random number k that represents a shared secret key for communication, the sending and receiving correspondents further agreeing on a base point on an elliptic polynomial (x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B))εEC^(xn+ny+2), wherein EC^(xn+ny+2) represents the elliptic polynomial, and further agreeing on a base point on the twist of the elliptic polynomial (x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB))εT EC^(xn+ny+2), wherein TEC^(nx+ny+2) represents the twist of the elliptic polynomial EC^(nx+ny+2), and wherein α is a non-quadratic residue element of F, the sending and receiving correspondents further agreeing on an initial vector S_(m) ⁻¹ , having a length of (nx+ny+1)N-bits;

the sending correspondent then performs the following steps:

c) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny) of a key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using a non-iterative embedding method to be described below;

d) if (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, a scalar multiplication is computed as (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=k(x _(0,TB) ,x _(1,TB) , . . . , x _(nx,TB) ,√{square root over ( α)}y _(0,B) ,y _(1,TB) , . . . , y _(ny,TB)) and (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,y _(0,k) ,y _(1,k) , . . . , y _(ny,k)), otherwise if α_(k)= α, then the scalar multiplication is computed as (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=k(x _(0,B) ,x _(1,B) , . . . , x _(nx,B) ,y _(0,B) ,y _(1,B) , . . . , y _(ny,B)) and (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,√{square root over ( α)}y _(0,k) ,y _(1,k) , . . . , y _(ny,k));

e) denoting the first (nx+ny+1)N-bit string of the message data to be s_(m) ₀ ′ and computing the (nx+ny+1)N-bit string s_(m) ₀ as s_(m) ₀ =s_(m) ₀ ′⊕s_(m) ⁻¹ ;

f) dividing the message bit string into M-bit strings, where ((nx+ny+1)N−L)>M>(N−L), wherein M and L are integers, and embedding the message bit string of the initial block of the message bit string s_(m) ₀ into the (nx+1) x-coordinates and ny y-coordinates of the elliptic message point (x_(0,m) ₀ , x_(1,m) ₀ , . . . , x_(nx,m) ₀ , √{square root over (α_(m) ₀ )}y_(0,m) ₀ , y_(1,m) ₀ , . . . , y_(ny,m) ₀ ) using the non-iterative embedding method;

g) if the message point of the 0-th block is on the elliptic polynomial, then the message authentication code points are computed as: (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x _(0,Tc) ₀ ,x _(1,Tc) ₀ , . . . , x _(nx,Tc) ₀ ,√{square root over ( α)}y _(0,Tc) ₀ ,y _(1,Tc) ₀ , . . . , y _(ny,Tc) ₀ )=(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )

Otherwise, they are computed as: (x _(0,Tc) ₀ ,x _(1,Tc) ₀ , . . . , x _(nx,Tc) ₀ ,√{square root over ( α)}y _(0,Tc) ₀ ,y _(1,Tc) ₀ , . . . , y _(ny,Tc) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,√{square root over ( α)}y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ ) and (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ );

h) initializing i as I=1, and repeating the following steps i) through l) until i>u:

i) denoting the (nx+ny+1)N-bit string of an i_(th) message data to be S_(m) _(i) and computing the (nx+ny+1)N-bit string s_(m) _(i) as s_(m) _(i) =s_(m) _(i) ′⊕s_(m) _(i−1) ;

j) embedding the message (nx+ny+1)N-bit string of an i-th block S_(m) _(i) ′ into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx) and the ny y-coordinates, y₁, . . . , y_(ny), of the elliptic message point (x_(0,m) _(i) , x_(1,m) _(i) , . . . , x_(nx,m) _(i) , √{square root over (α_(m) _(i) )}y_(0,m) _(i) , y_(1,m) _(i) , . . . , y_(ny,m) _(i) ) using the non-iterative embedding method;

k) doubling the points (x_(0,c) _(i−1) , x_(1,c) _(i−1) , . . . , x_(nx,c) _(i−1) , y_(0,c) _(i−1) , y_(1,c) _(i−1) , . . . , y_(ny,c) _(i−1) ) and (x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , y_(ny,Tc) _(i−1) ) as: (x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , y _(ny,c) _(i−1) )=2(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , y _(ny,c) _(i−1) ) and (x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , y _(ny,Tc) _(i−1) )=2(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , y _(ny,Tc) _(i−1) );

l) if the message point of the i-th block is on the elliptic polynomial, then the message authentication code points are computed as: (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i) , . . . , x _(ny,c) _(i−1) ) and (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,√{square root over ( α)}y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , x _(ny,Tc) _(i−1) ).

Otherwise, they are computed as: (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,√{square root over ( α)}y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,√{square root over ( α)}y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , x _(ny,Tc) _(i−1) ) and (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i) , . . . , x _(ny,c) _(i−1) );

m) a set of appropriate bits of the x-coordinates x_(0,Tc), x_(1,Tc), . . . , x_(nx,Tc), y-coordinates y_(1,c), . . . , y_(ny,c) and a sign bit of y_(0,c) of the message authentication code point (x_(0,c) _(u) , x_(1,c) _(u) , . . . , x_(nx,c) _(u) , y_(0,c) _(u) , y_(1,c) _(u) , . . . , y_(ny,c) _(u) ), and a set of appropriate bits of the x-coordinates x_(0,Tc), x_(1,Tc), . . . , x_(nx,Tc), y-coordinates y_(1,Tc), . . . , y_(ny,Tc) and the sign bit of y_(0,Tc) of the message authentication code point (x_(0,Tc) _(u) , x_(1,Tc) _(u) , . . . , x_(nx,Tc) _(u) , √{square root over ( α)}y_(0,Tc) _(u) , y_(1,Tc) _(u) , . . . , y_(ny,Tc) _(u) ) are concatenated together to form the message authentication code, which is then appended to the message to be sent together with, if necessary, any additional information needed to help authenticate the message bit string at the receiving correspondent without compromising security and sent to the receiving correspondent;

the receiving correspondent then performs the following steps:

n) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny) of the key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using the non-iterative embedding method;

o) if the point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, then a scalar multiplication is performed as (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=k(x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB)) and (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), y_(0,k), y_(1,k), . . . , y_(ny,k)), otherwise if α_(k)=α_(o) then the scalar multiplication is performed as (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=k(x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B)) (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), y_(0,k), y_(1,k), . . . , y_(ny,k));

p) denoting the first (nx+ny+1)N-bit string of the message data to be s_(rm) ₀ ′ and computing the (nx+ny+1)N-bit string s_(rm) ₀ ) as s_(rm) ₀ =s_(rm) ₀ ′⊕s_(rm) ⁻¹ ;

q) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), wherein M and L are integers, and embedding the message bit string of the initial block of the message bit string s_(rm) ₀ into the (nx+1) x-coordinates and the ny y-coordinates of the received message elliptic point (x_(0,rm) ₀ , x_(1,rm) ₀ , . . . , x_(nx,rm) ₀ , √{square root over (α_(rm) ₀ )}y_(0,rm) ₀ , y_(1,rm) ₀ , . . . , y_(ny,rm) ₀ ) using the non-iterative embedding method;

r) if the message point of the 0-th block is on the elliptic polynomial, then the message authentication points are computed as: (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x _(0,rTc) ₀ ,x _(1,rTc) ₀ , . . . , x _(nx,rTc) ₀ ,√{square root over ( α)}y _(0,rTc) ₀ ,y _(1,rTc) ₀ , . . . , y _(ny,rTc) ₀ )=(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ ).

Otherwise they are computed as: (x _(0,rTc) ₀ ,x _(1,rTc) ₀ , . . . , x _(nx,rTc) ₀ ,√{square root over ( α)}y _(0,rTc) ₀ ,y _(1,rTc) ₀ , . . . , y _(ny,rTc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,rTS) ₀ ,x _(1,rTS) ₀ , . . . , x _(nx,rTS) ₀ ,√{square root over ( α)}y _(0,rTS) ₀ ,y _(1,rTS) ₀ , . . . , x _(ny,rTS) ₀ ) and (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )

s) initializing the integer i as i=0, the following steps t) through w) are repeated until i>u:

t) denoting the (nx+ny+1)N-bit string of an i_(th) message data to be s_(rm) _(i) ′ and computing the (nx+ny+1)N-bit string s_(rm) _(i) as s_(rm) _(i) =s_(rm) _(i) ′⊕s_(rm) _(i−1) ;

u) embedding the received message (nx+ny+1)N-bit string of the i-th block s_(rm) _(i) the into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of the received message elliptic point (x_(0,rm) _(i) , x_(1,rm) _(i) , . . . , x_(nx,rm) _(i) , √{square root over (α_(rm) _(i) )}y_(0,rm) _(i) , y_(1,rm) _(i) , . . . , y_(ny,rm) _(i) ) using the non-iterative embedding method;

v) doubling the points (x_(0,rc) _(i−1) , x_(1,rc) _(i−1) , . . . , x_(nx,rc) _(i−1) , y_(0,rc) _(i−1) , y_(1,rc) _(i−1) , . . . , y_(ny,rc) _(i−1) ) and (x_(0,rTc) _(i−1) , x_(1,rTc) _(i−1) , . . . , x_(nx,rTc) _(i−1) , y_(0,rTc) _(i−1) , y_(1,rTc) _(i−1) , . . . , y_(ny,rTc) _(i−1) ) as: (x _(0,rc) _(i−1) , x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) , y _(0,rc) _(i−1) , y _(1,rc) _(i−1) , . . . , y _(ny,rc) _(i−1) )=2(x _(0,rc) _(i−1) , x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) , y _(0,rc) _(i−1) , y _(1,rc) _(i−1) , . . . , y _(ny,rc) _(i−1) ) (x _(0,rTc) _(i−1) , x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) , y _(0,rTc) _(i−1) , y _(1,rTc) _(i−1) , . . . , y _(ny,rTc) _(i−1) )=2(x _(0,rTc) _(i−1) , x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) , y _(0,rTc) _(i−1) , y _(1,rTc) _(i−1) , . . . , y _(ny,rTc) _(i−1) );

w) if the received message point of the i-th block is on the elliptic polynomial, then the message authentication points are computed as (x _(0,rc) _(i) , x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) , y _(0,rc) _(i) , y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rm) _(i) , x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) , y _(0,rm) _(i) , y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rc) _(i−1) , x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) , y _(0,rc) _(i−1) , y _(1,rc) _(i) , . . . , x _(ny,rc) _(i−1) ) and (x_(0,rTc) _(i) , x_(1,rTc) _(i) , . . . , x_(nx,rTc) _(i) , √{square root over ( α)}y_(0,rTc) _(i) , y_(1,rTc) _(i) , . . . , y_(ny,rTc) _(i) )=(x_(0,rTc) _(i−1) , x_(1,rTc) _(i−1) , . . . , x_(nx,rTc) _(i−1) , √{square root over ( α)}y_(0,rTc) _(i−1) , y_(1,rTc) _(i−1) , . . . , x_(ny,rTc) _(i−1) ). Otherwise they are computed as: (x _(0,rTc) _(i) ,x _(1,rTc) _(i) , . . . , x _(nx,rTc) _(i) ,√{square root over ( α)}y _(0,rTc) _(i) ,y _(1,rTc) _(i) , . . . , y _(ny,rTc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,√{square root over ( α)}y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rTc) _(i−1) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) ,√{square root over ( α)}y _(0,rTc) _(i−1) ,y _(1,rTc) _(i−1) , . . . , x _(ny,rTc) _(i−1) ) and (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i) , . . . , x _(ny,rc) _(i−1) );

x) the appropriate bits of the x-coordinates x_(0,rc), x_(1,rc), . . . , x_(nx,rc), y-coordinates y_(1,rc), . . . , y_(ny,rc) and the sign bit of y_(0,rc) of the message authentication code point (x_(0,rc) _(u) , x_(1,rc) _(u) , . . . , x_(nx,rc) _(u) , y_(0,rc) _(u) , y_(1,rc) _(u) , . . . , y_(ny,rc) _(u) ), and the appropriate bits of the x-coordinates x_(0,rTc), x_(1,rTc), . . . , x_(nx,rTc), y-coordinates y_(1,rTc), . . . , y_(ny,rTc) and the sign bit of y_(0,rTc) of the message authentication point (x_(0,rTc) _(u) , x_(1,rTc) _(u) , . . . , x_(nx,rTc) _(u) , √{square root over ( α)}y_(0,rTc) _(u) , y_(1,rTc) _(u) , . . . , y_(ny,rTc) _(u) ) are concatenated together to form the message authentication code of the received message; and

y) if the received message authentication code is equal to the message authentication code of the received message generated at the receiving correspondent, the received message is authenticated.

The above methods can also be applied to find the MAC for bit streams of media data such as text, audio, video, or a multimedia data. Further, the MAC methods given above can be applied to the verification of media data. For efficiency of computation, the above methods are used with a pre-processing stage, which is used to compress the media data prior to the application of the MAC method. Either a lossless compression method or a lossy compression method can be used to compress the media data in the pre-processing stage. The bit string of the compressed message at the output of the pre-processing stage is then used as the input to the MAC generation methods.

With regard to data embedding, in the following, it is assumed that the maximum block size that can be embedded into the (nx+1) x-coordinates and the ny y-coordinates is (nx+ny+1)N bits, and that the compressed media data bit string length is a multiple of (nx+ny+1)N, such as (u+1) (nx+ny+1)N. In other words, the number of (nx+ny+1)N-bit blocks in a message bit string is (u+1). The data embedding is then performed as follows:

a) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), wherein M and L are integers;

b) dividing each M-bit string into (nx+ny+1) strings mx₀, mx₁, . . . , mx_(nx), my₁, . . . , my_(ny), wherein the length of string mx₀ is less than or equal to (N−L) bits, the length of each of the remaining (nx+ny) strings mx₁, . . . , mx_(nx), my₁, . . . , my_(ny) being less than or equal to N bits, wherein each of the (nx+ny+1) bit strings has an equivalent value which is an element of the finite field F;

c) assigning the value of the bit strings of mx₁, . . . , mx_(nx) to x₁, . . . , x_(nx);

d) assigning the value of the bit strings of my₁, . . . , my_(ny) to y₁, . . . , y_(ny);

e) assigning the value of the bit string of mx₀ to x₀;

f) substituting the values of x₀, x₁, . . . , x_(nx) and y₁, . . . , y_(ny) in a selected elliptic polynomial equation with (nx+1) x-coordinates and (ny+1) y-coordinates to form a quadratic equation in y₀ of the form y₀ ²+Ay₀+B=0, wherein A and B are variable coefficients;

g) if the quadratic equation in y_(o) has solutions y _(0,1), y _(0,2), wherein the solutions are elements of the finite field F, then a method agreed upon by the sending and receiving correspondent is used to select one of the solutions y _(0,1), y _(0,2), and further assigning the chosen solution to y_(o), assigning α_(m)=1, and the message point is given by (x_(0,m), x_(1,m), . . . , x_(nx,m), y_(0,m), y_(1,m), . . . , y_(ny,m)), otherwise a method agreed upon by the sending and receiving correspondent is used to select one of the solutions √{square root over ( α)} y _(0,1), √{square root over ( α)} y _(0,2), and then further assign the chosen solution to y_(o) assigning α_(m)= α and the message point is given by (x_(0,m), x_(1,m), . . . , x_(nx,m), √{square root over ( α)}y_(0,m), y_(1,m), . . . , y_(ny,m))

h) denoting the message point as: (x_(0,m), x_(1,m), . . . , x_(nx,m), √{square root over (α)}y_(0,m), y_(1,m), . . . , y_(ny,m)), where the point is on the elliptic polynomial if α_(m) _(i) =1, and the point is on the twist if α_(m) _(i) = α.

It should be noted that in the case of quadratic equations with A=0, the Legendre symbol can be used as an efficient test for the existence of a solution of the quadratic equation above. The Legendre Symbol is used to test whether an element of F(p) has a square root or not, i.e., whether an element is quadratic residue or not. The Legendre Symbol and test are as follows. Given an element of a finite field F(p), such as d, the Legendre symbol is defined as

$\left( \frac{d}{p} \right).$ In order to test whether d is quadratic residue or not, the Legendre symbol,

$\left( \frac{d}{p} \right),$ is computed such that:

$\left( \frac{d}{p} \right) = \left\{ {\begin{matrix} {+ 1} & {{if}\mspace{14mu} x\mspace{14mu}{is}\mspace{14mu}{quadratic}\mspace{14mu}{residue}} \\ 0 & {{{if}\mspace{14mu} x} \equiv {0\mspace{14mu}{mod}\mspace{14mu}{F(p)}}} \\ {- 1} & {otherwise} \end{matrix}.} \right.$

In the above, the MACs use the scalar multiplication k_(m)(x_(Pu),y_(Pu)). It should be noted that, in order to find a collision means, that there are two message bits strings m and m′ such that their integer values k_(m) and k_(m′) will lead to k_(m)k(x_(B),y_(B))≡k_(m′)k(x_(B),y_(B)). This collision implies that integers can be found such that k_(m)k−k_(m′)k=l*#EC, where #EC is the order of the group (EC,+). This is equivalent to solving the elliptic polynomial discrete logarithm problem. This also applies to finding a collision for the points on the twist of an elliptic polynomial, k_(m)k(x_(TB), √{square root over ( α)}y_(TB))=k_(m′)k(x_(TB), √{square root over ( α)}y_(TB)).

Thus, security of the MACs depends on the security of the underlying elliptic polynomial cryptography. The security of elliptic polynomial cryptosystems is assessed by both the effect on the solution of the elliptic polynomial discrete logarithmic problem (ECDLP) and power analysis attacks.

It is well known that the elliptic polynomial discrete logarithm problem (ECDLP) is apparently intractable for non-singular elliptic polynomials. The ECDLP problem can be stated as follows: given an elliptic polynomial defined over F that needs N-bits for the representation of its elements, an elliptic polynomial point (x_(p),y_(p))εEC, defined in affine coordinates, and a point (x_(Q),y_(Q))εEC, defined in affine coordinates, determine the integer k, 0≦k≦#F, such that (x_(Q),y_(Q))=k(x_(p),y_(p)), provided that such an integer exists. In the below, it is assumed that such an integer exists.

The most well known attack used against the ECDLP is the Pollard ρ-method, which has a complexity of O(√{square root over (πK)}/2), where K is the order of the underlying group, and the complexity is measured in terms of an elliptic polynomial point addition.

Since the underlying cryptographic problems used in the above block cipher chaining methods is the discrete logarithm problem, which is a known difficult mathematical problem, it is expected that the security of the above methods are more secure than prior art ciphers which are not based on such a mathematically difficult problem.

Projective coordinate can also be used by the sending correspondent and the receiving correspondent to embed extra message data bits in the projective coordinate wherein the addition of the corresponding elliptic points is defined in (nx+ny+3) dimensional space where there are (nx+1) x-coordinates, (ny+1) y-coordinates and one projective coordinate.

The equations for the addition rule can be obtained by using the elliptic polynomial equation with (nx+1) x-coordinates and (nx+1) y coordinates in projective coordinates and substituting a straight line equation to obtain a cubic equation in terms of one of the x-coordinates. This cubic equation can be used to identify the third point of intersection between a straight line and the elliptic polynomial in (nx+ny+3) dimensions given two other intersection points. This third point of intersection is used to identify the sum of the given two points.

It will be understood that the MACs based on elliptic polynomial cryptography described above may be implemented by software stored on a medium readable by a computer and executing as set of instructions on a processor (including a microprocessor, microcontroller, or the like) when loaded into main memory in order to carry out a cryptographic system of secure communications in a computer network. As used herein, a medium readable by a computer includes any form of magnetic, optical, mechanical, laser, or other media readable by a computer, including floppy disks, hard disks, compact disks (CDs), digital versatile disk (DVD), laser disk, magnetic tape, paper tape, punch cards, flash memory, etc.

It is to be understood that the present invention is not limited to the embodiments described above, but encompasses any and all embodiments within the scope of the following claims. 

We claim:
 1. A computerized method of generating an elliptic polynomial-based message authentication code, comprising the steps of: a) defining a maximum block size that can be embedded into (nx+1) x-coordinates and ny y-coordinates, wherein n is an integer, and setting the maximum block size to be (nx+ny+1)N bits, wherein N is an integer; b) a sending correspondent and a receiving correspondent agree upon the values of nx and ny, and further agree on a set of coefficients a_(1k), a_(2kl), a_(3k), c_(1lki), c_(2kl), c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k) & b_(c)εF, wherein F represents a finite field where the field's elements can be represented in N-bits, the sending and receiving correspondents further agreeing on a random number k that represents a shared secret key for communication, the sending and receiving correspondents further agreeing on a base point on an elliptic polynomial (x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B))εEC^(xn+ny+2), wherein EC^(nx+ny+2) represents the elliptic polynomial, and further agreeing on a base point on the twist of the elliptic polynomial (x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB))εEC^(xn+ny+2), wherein TEC^(nx+ny+2) represents the twist of the elliptic polynomial EC^(nx+ny+2), and further, α is a non-quadratic residue element of F; the sending correspondent then performs the following steps: c) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of a key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using a non-iterative embedding method; d) if (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, a scalar multiplication is computed as: (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=k(x _(0,TB) ,x _(1,) , . . . , x _(nx,TB) ,√{square root over ( α)}y _(0,TB) ,y _(1,TB) , . . . , y _(ny,TB)) and (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,y _(0,k) ,y _(1,k) , . . . , y _(ny,k)), otherwise, if α_(k)= α, then the scalar multiplication is computed as: (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=k(x _(0,B) ,x _(1,B) , . . . , x _(nx,B) ,y _(0,B) ,y _(1,B) , . . . , y _(ny,B)) and (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,√{square root over ( α)}y _(0,k) ,y _(1,k) , . . . , y _(ny,k)); e) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), where M and L are integers, and embedding the message bit string of the initial block of the message bit string into the (nx+1) x-coordinates and the ny y-coordinates of the elliptic message point (x_(0,m) ₀ , x_(1,m) ₀ , . . . , x_(nx,m) ₀ , √{square root over (α_(m) ₀ )}y_(0,m) ₀ , y_(1,m) ₀ , . . . , y_(ny,m) ₀ ) using the non-iterative embedding method; f) if the message point of the 0-th block is on the elliptic polynomial, then the message authentication code points are computed as (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,√{square root over ( α)}y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,√{square root over ( α)}y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,√{square root over ( α)}y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x_(0,Tc) ₀ , x_(1,Tc) ₀ , . . . , x_(nx,Tc) ₀ , √{square root over ( α)}y_(0,Tc) ₀ , y_(1,Tc) ₀ , . . . , y_(ny,Tc) ₀ )=(x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ ), otherwise, they are computed as: (x _(0,Tc) ₀ ,x _(1,Tc) ₀ , . . . , x _(nx,Tc) ₀ ,y _(0,Tc) ₀ ,y _(1,Tc) ₀ , . . . , y _(ny,Tc) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ ) and (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ); g) initializing an integer i as i=1, and repeating the following steps h) and i) and incrementing i at each step until all of the message data blocks are processed: h) embedding the message (nx+ny+1)N-bit string of an i-th block into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx) and the ny y-coordinates, y₁, . . . , y_(ny), of the elliptic message point (x_(0,m) _(i) , x_(1,m) _(i) , . . . , x_(nx,m) _(i) , √{square root over (α_(m) _(i) )}y_(0,m) _(i) , y_(1,m) _(i) , . . . , y_(ny,m) _(i) ) using the non-iterative embedding method; i) if the message point of the i-th block is on the elliptic polynomial, then the message authentication code points are computed as (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i) , . . . , x _(ny,c) _(i−1) ) and (x_(0,Tc) _(i) , x_(1,Tc) _(i) , . . . , x_(nx,Tc) _(i) , √{square root over ( α)}y_(0,Tc) _(i) , y_(1,Tc) _(i) , . . . , y_(ny,Tc) _(i) )=(x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , √{square root over ( α)}y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , x_(ny,Tc) _(i−1) ), otherwise, they are computed as: (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,√{square root over ( α)}y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,√{square root over ( α)}y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , x _(ny,Tc) _(i−1) ) and (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i) , . . . , x _(ny,c) _(i−1) ); j) a set of appropriate bits of the x-coordinates x_(0,c), x_(1,c), . . . , x_(nx,c), y-coordinates y_(1,c), . . . , y_(ny,c) and a sign bit of y_(0,c) of the message authentication code point (x_(0,c) _(u) , x_(1,c) _(u) , . . . , x_(nx,c) _(u) , y_(0,c) _(u) , y_(1,c) _(u) , . . . , y_(ny,c) _(u) ) and a set of appropriate bits of the x-coordinates x_(0,Tc), x_(1,Tc), . . . , x_(nx,Tc), y-coordinates y_(1,Tc), . . . , y_(ny,Tc) and the sign bit of y_(0,Tc) of the message authentication code point (x_(0,Tc) _(u) , x_(1,Tc) _(u) , . . . , x_(nx,Tc) _(u) , √{square root over ( α)}y_(0,Tc) _(u) , y_(1,Tc) _(u) , . . . , y_(ny,Tc) _(u) ) are concatenated together to form the message authentication code, which is then appended to the message to be sent together with, if necessary, any additional information needed to help authenticate the message bit string at the receiving correspondent without compromising security and sent to the receiving correspondent; the receiving correspondent then performs the following steps: k) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of the key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using the non-iterative embedding method; l) if the point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, then a scalar multiplication is performed as (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=k(x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,TB), y_(1,TB), . . . , y_(ny,TB)) and (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), y_(0,k), y_(1,k), . . . , y_(ny,k)) otherwise i if α_(k)=α_(o) then the scalar multiplication is performed as (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=k(x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B)) and (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)); m) dividing the received message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), where M and L are integers, and embedding the received message (nx+ny+1)N-bit string of the initial block of the received message bit string into the (nx+1) x-coordinates and the ny y-coordinates of the received message elliptic point (x_(0,rm) ₀ , x_(1,rm) ₀ , . . . , x_(nx,rm) ₀ , √{square root over (α_(rm) ₀ )}y_(0,rm) ₀ , y_(1,rm) ₀ , . . . , y_(ny,rm) ₀ ) using the non-iterative embedding method; n) if the message point of the initial block is on the elliptic polynomial, then the message authentication points are computed as: (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x_(0,rTc) ₀ , x_(1,rTc) ₀ , . . . , x_(nx,rTc) ₀ , √{square root over ( α)}y_(0,rTc) ₀ , y_(1,rTc) ₀ , . . . , y_(ny,rTc) ₀ )=(x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ ), otherwise, they are computed as: (x _(0,rTc) ₀ ,x _(1,rTc) ₀ , . . . , x _(nx,rTc) ₀ ,√{square root over ( α)}y _(0,rTc) ₀ ,y _(1,rTc) ₀ , . . . , y _(ny,rTc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,√{square root over ( α)}y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,rTS) ₀ ,x _(1,rTS) ₀ , . . . , x _(nx,rTS) ₀ ,√{square root over ( α)}y _(0,rTS) ₀ ,y _(1,rTS) ₀ , . . . , x _(ny,rTS) ₀ ) and (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ); o) initializing the integer i as i=0, and repeating the following steps p) and q) and incrementing i at each step until all of the message data blocks are processed: p) embedding the received message (nx+ny+1)N-bit string of the i-th block the into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny)) of the received message elliptic point (x_(0,rm) _(i) , x_(1,rm) _(i) , . . . , x_(nx,rm) _(i) , √{square root over (α_(rm) _(i) )}y_(0,rm) _(i) , y_(1,rm) _(i) , . . . , y_(ny,rm) _(i) ) using the non-iterative embedding method; q) if the received message point of the i-th block is on the elliptic polynomial, then the message authentication points are computed as: (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i−1) , . . . , x _(ny,rc) _(i−1) ) and (x_(0,rTc) _(i) , x_(1,rTc) _(i) , . . . , x_(nx,rTc) _(i) , √{square root over ( α)}y_(0,rTc) _(i) , y_(1,rTc) _(i) , . . . , y_(ny,rTc) _(i) )=(x_(0,rTc) _(i−1) , x_(1,rTc) _(i−1) , . . . , x_(nx,rTc) _(i−1) , √{square root over ( α)}y_(0,rTc) _(i−1) , y_(1,rTc) _(i−1) , . . . , x_(ny,rTc) _(i−1) ), otherwise, they are computed as: (x _(0,rTc) _(i) ,x _(1,rTc) _(i) , . . . , x _(nx,rTc) _(i) ,√{square root over ( α)}y _(0,rTc) _(i) ,y _(1,rTc) _(i) , . . . , y _(ny,rTc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,√{square root over ( α)}y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rTc) _(i−1) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) ,√{square root over ( α)}y _(0,rTc) _(i−1) ,y _(1,rTc) _(i) , . . . , x _(ny,rTc) _(i−1) ); r) the appropriate bits of the x-coordinates x_(0,rc), x_(1,rc), . . . , x_(nx,rc), y-coordinates y_(1,rc), . . . , y_(ny,rc) and the sign bit of y_(0,rc) of the message authentication code point (x_(0,rc) _(u) , x_(1,rc) _(u) , . . . , x_(nx,rc) _(u) , y_(0,rc) _(u) , y_(1,rc) _(u) , . . . , y_(ny,rc) _(u) ), and the appropriate bits of the x-coordinates x_(0,rTc), x_(1,rTc), . . . , x_(nx,rTc), y-coordinates y_(1,rTc), . . . , y_(ny,rTc) and the sign bit of y_(0,rTc) of the message authentication point (x_(0,rTc) _(u) , x_(1,rTc) _(u) , . . . , x_(nx,rTc) _(u) , y_(0,rTc) _(u) , y_(1,rTc) _(u) , . . . , y_(ny,rTc) _(u) ) are concatenated together to form the message authentication code of the received message; and s) if the received message authentication code is equal to the message authentication code of the received message generated at the receiving correspondent, the received message is authenticated.
 2. The computerized method of generating an elliptic polynomial-based message authentication code as recited in claim 1, wherein the non-iterative embedding method comprises the steps of: a) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), wherein M and L are integers; b) dividing each M-bit string into (nx+ny+1) strings mx₀, mx₁, . . . , mx_(nx), my₁, . . . , my_(ny), wherein the length of string mx₀ is less than or equal to (N−L) bits, the length of each of the remaining (nx+ny) strings mx₁, . . . , mx_(nx), my₁, . . . , my_(ny) being less than or equal to N bits, wherein each of the (nx+ny+1) bit strings has an equivalent value which is an element of the finite field F; c) assigning the value of the bit strings of mx₁, . . . , mx_(nx) to x₁, . . . , x_(nx); d) assigning the value of the bit strings of my₁, . . . , my_(ny) to y₁, . . . , y_(ny); e) assigning the value of the bit string of mx₀ to x₀; f) substituting the values of x₀, x₁, . . . , x_(nx) and y₁, . . . , y_(ny) in a selected elliptic polynomial equation with (nx+1) x-coordinates and (ny+1) y-coordinates to form a quadratic equation in y₀ of the form y₀ ²+Ay₀+B=0, wherein A and B are variable coefficients; g) if the quadratic equation in y_(o) has solutions y _(0,1), y _(0,2), wherein the solutions are elements of the finite field F, then a method agreed upon by the sending and receiving correspondent is used to select one of the solutions y _(0,1), y _(0,2), and further assigning the chosen solution to y_(o), assigning α_(m)=1, and the message point is given by (x_(0,m), x_(1,m), . . . , x_(nx,m), y_(0,m), y_(1,m), . . . , y_(ny,m)), otherwise a method agreed upon by the sending and receiving correspondent is used to select one of the solutions √{square root over ( α)} y _(0,1), √{square root over ( α)} y _(0,2), and then further assign the chosen solution to y_(o), assigning α_(m)= α and the message point is given by (x_(0,m), x_(1,m), . . . , x_(nx,m), √{square root over ( α)}y_(0,m), y_(1,m), . . . , y_(ny,m)); and h) denoting the message point as (x_(0,m), x_(1,m), . . . , x_(nx,m), √{square root over (α)}y_(0,m), y_(1,m), . . . , y_(ny,m)).
 3. A computerized method of generating an elliptic polynomial-based message authentication code, comprising the steps of: a) defining a maximum block size that can be embedded into (nx+1) x-coordinates and ny y-coordinates, wherein n is an integer, and setting the maximum block size to be (nx+ny+1)N bits, wherein N is an integer; b) a sending correspondent and a receiving correspondent agree upon the values of nx and ny, where nx and ny are integers, and further agree on a set of coefficients a_(1k), a_(2kl), a_(3k), c_(1lki), c_(2kl), c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k) & b_(c)εF, wherein F represents a finite field and where the elements of F are represented in N-bits, the sending and receiving correspondents further agreeing on a random number k that represents a shared secret key for communication, the sending and receiving correspondents further agreeing on a base point on an elliptic polynomial (x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B))εEC^(xn+ny+2), wherein EC^(nx+ny+2) represents the elliptic polynomial, and further agreeing on a base point on the twist of the elliptic polynomial (x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB))εTEC^(xn+ny+2), wherein TEC^(nx+ny+2) represents the twist of the elliptic polynomial EC^(nx+ny+2), and where α is a non-quadratic residue element of F; the sending correspondent then performs the following steps: c) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny)) of a key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using a non-iterative embedding method; d) if (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, a scalar multiplication is computed as: (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=k(x _(0,TB) ₀ ,x _(1,TB) ₀ , . . . , x _(nx,TB) ₀ ,√{square root over ( α)}y _(0,B) ₀ ,y _(1,TB) ₀ , . . . , y _(ny,TB) ₀ ) and (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), y_(0,k), y_(1,k), . . . , y_(ny,k)), otherwise, if α_(k)= α, then the scalar multiplication is computed as: (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=k(x _(0,B) ,x _(1,B) , . . . , x _(nx,B) ,y _(0,B) ,y _(1,B) , . . . , y _(ny,B)) and (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,y _(0,k) ,y _(1,k) , . . . , y _(ny,k)); e) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), where M and L are integers, and embedding the message bit string of the initial, block of the message bit string into the (nx+1) x-coordinates and the ny y-coordinates of the elliptic message point (x_(0,m) ₀ , x_(1,m) ₀ , . . . , x_(nx,m) ₀ , √{square root over (α_(m) ₀ )}y_(0,m) ₀ , y_(1,m) ₀ , . . . , y_(ny,m) ₀ ) using the non-iterative embedding method; f) if the message point of the 0-th block is on the elliptic polynomial, then the message authentication code points are computed as: (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , y _(ny,S) ₀ ) and (x_(0,Tc) ₀ , x_(1,Tc) ₀ , . . . , x_(nx,Tc) ₀ , √{square root over ( α)}y_(0,Tc) ₀ , y_(1,Tc) ₀ , . . . , y_(ny,Tc) ₀ )=(x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , y_(ny,TS) ₀ ), otherwise, they are computed as: (x _(0,Tc) ₀ ,x _(1,Tc) ₀ , . . . , x _(nx,Tc) ₀ ,√{square root over ( α)}y _(0,Tc) ₀ ,y _(1,Tc) ₀ , . . . , y _(ny,Tc) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,√{square root over ( α)}y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , y _(ny,TS) ₀ ) and (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ); g) initializing an integer i as i=1, and repeating the following steps h) through j) and incrementing i at each step until all of the message data blocks are processed: h) embedding the message (nx+ny+1)N-bit string of an i-th block into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx) and the ny y-coordinates, y₁, . . . , y_(ny), of the elliptic message point (x_(0,m) _(i) , x_(1,m) _(i) , . . . , x_(nx,m) _(i) , √{square root over (α_(m) _(i) )}y_(0,m) _(i) , y_(1,m) _(i) , . . . , y_(ny,m) _(i) ) using the non-iterative embedding method; i) doubling the points (x_(0,c) _(i−1) , x_(1,c) _(i−1) , . . . , x_(nx,c) _(i−1) , y_(0,c) _(i−1) , y_(1,c) _(i−1) , . . . , y_(ny,c) _(i−1) ) and (x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , y_(ny,Tc) _(i−1) ) as: (x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , y _(ny,c) _(i−1) )=2(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , y _(ny,c) _(i−1) ) and (x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , y _(ny,Tc) _(i−1) )=2(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , y _(ny,Tc) _(i−1) ); j) if the message point of the i-th block is on the elliptic polynomial, then the message authentication code points are computed as: (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , x _(ny,c) _(i−1) ) and (x_(0,Tc) _(i) , x_(1,Tc) _(i) , . . . , x_(nx,Tc) _(i) , √{square root over ( α)}y_(0,Tc) _(i) , y_(1,Tc) _(i) , . . . , y_(ny,Tc) _(i) )=(x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , √{square root over ( α)}y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , x_(ny,Tc) _(i−1) ), otherwise, they are computed as: (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,√{square root over ( α)}y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,√{square root over ( α)}y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , x _(ny,Tc) _(i−1) ) and (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , x _(ny,c) _(i−1) ); k) a set of appropriate bits of the x-coordinates x_(0,c), x_(1,c), . . . , x_(nx,c), y-coordinates y_(1,c), . . . , y_(ny,c) and a sign bit of y₀,c of the message authentication code point (x_(0,c) _(u) , x_(1,c) _(u) , . . . , x_(nx,c) _(u) , y_(0,c) _(u) , y_(1,c) _(u) , . . . , y_(ny,c) _(u) ) and a set of appropriate bits of the x-coordinates x_(0,Tc), x_(1,Tc), . . . , x_(nx,Tc), y-coordinates y_(1,Tc), . . . , y_(ny,Tc) and the sign bit of y_(0,Tc) of the message authentication code point (x_(0,Tc) _(u) , x_(1,Tc) _(u) , . . . , x_(nx,Tc) _(u) , √{square root over ( α)}y_(0,Tc) _(u) , y_(1,Tc) _(u) , . . . , y_(ny,Tc) _(u) ) are concatenated together to form the message authentication code, which is then appended to the message to be sent together with, if necessary, any additional information needed to help authenticate the message bit string at the receiving correspondent without compromising security and sent to the receiving correspondent; the receiving correspondent then performs the following steps: l) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of the key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using the non-iterative embedding method; m) if the point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, then a scalar multiplication is performed as (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=k(x_(0,TB) ₀ , x_(1,TB) ₀ , . . . , x_(nx,TB) ₀ , √{square root over ( α)}y_(0,B) ₀ , y_(1,TB) ₀ , . . . , y_(ny,TB) ₀ ) and (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), y_(0,k), y_(1,k), . . . , y_(ny,k)) otherwise if α_(k)= α then the scalar multiplication is performed as (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=k(x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B)) and (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over ( α)}y_(0,k), y_(1,k), . . . , y_(ny,k)); n) dividing the received message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), where M and L are integers, and embedding the received message bit string of the initial block of the received message bit string into the (nx+1) x-coordinates and the ny y-coordinates of the go received message elliptic point (x_(0,rm) ₀ , x_(1,rm) ₀ , . . . , x_(nx,rm) ₀ , √{square root over (α_(rm) ₀ )}y_(0,rm) ₀ , y_(1,rm) ₀ , . . . , y_(ny,rm) ₀ ) using the non-iterative embedding method; o) if the message point of the initial-block is on the elliptic polynomial, then the message authentication points are computed as: (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x _(0,rTc) ₀ ,x _(1,rTc) ₀ , . . . , x _(nx,rTc) ₀ , √{square root over ( α)}y _(0,rTc) ₀ ,y _(1,rTc) ₀ , . . . , y _(ny,rTc) ₀ )=(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ , √{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ ), otherwise they are computed as: (x _(0,rTc) ₀ ,x _(1,rTc) ₀ , . . . , x _(nx,rTc) ₀ ,√{square root over ( α)}y _(0,rTc) ₀ ,y _(1,rTc) ₀ , . . . , y _(ny,rTc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,√{square root over ( α)}y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,rTS) ₀ ,x _(1,rTS) ₀ , . . . , x _(nx,rTS) ₀ ,√{square root over ( α)}y _(0,rTS) ₀ ,y _(1,rTS) ₀ , . . . , x _(ny,rTS) ₀ ) and (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ); p) initializing the integer i as i=0, and repeating the following steps q) through s) and incrementing i at each step until all of the message data blocks are processed: q) embedding the received message (nx+ny+1)N-bit string of the i-th block the into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of the received message elliptic point (x_(0,rm) _(i) , x_(1,rm) _(i) , . . . , x_(nx,rm) _(i) , √{square root over (α_(rm) _(i) )}y_(0,rm) _(i) , y_(1,rm) _(i) , . . . , y_(ny,rm) _(i) ) using the non-iterative embedding method; r) doubling the points (x_(0,S) _(i) , x_(1,S) _(i) , . . . , x_(nx,S) _(i) , y_(0,S) _(i) , y_(1,S) _(i) , . . . , x_(ny,S) _(i) ) and (x_(0,TS) _(i) , x_(1,TS) _(i) , . . . , x_(nx,TS) _(i) , √{square root over ( α)}y_(0,TS) _(i) ,y_(1,TS) _(i) , . . . , x_(ny,TS) _(i) ) as: (x _(0,S) _(i) ,x _(1,S) _(i) , . . . , x _(nx,S) _(i) ,y _(0,S) _(i) ,y _(1,S) _(i) , . . . , x _(ny,S) _(i) )=2(x _(0,S) _(i−1) ,x _(1,S) _(i−1) , . . . , x _(nx,S) _(i−1) ,y _(0,S) _(i−1) ,y _(1,S) _(i−1) , . . . , x _(ny,S) _(i−1) ) (x _(0,TS) _(i) ,x _(1,TS) _(i) , . . . , x _(nx,TS) _(i) ,√{square root over ( α)}y _(0,TS) _(i) ,y _(1,TS) _(i) , . . . , x _(ny,TS) _(i) )=2(x _(0,TS) _(i−1) ,x _(1,TS) _(i−1) , . . . , x _(nx,TS) _(i−1) ,√{square root over ( α)}y _(0,TS) _(i−1) ,y _(1,TS) _(i−1) , . . . , x _(ny,TS) _(i−1) ) s) if the received message point of the i-th block is on the elliptic polynomial, then the message authentication points are computed as (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i) , . . . , x _(ny,rc) _(i−1) ) and (x _(0,rTc) _(i) ,x _(1,rTc) _(i) , . . . , x _(nx,rTc) _(i) , √{square root over ( α)}y _(0,rTc) _(i) ,y _(1,rTc) _(i) , . . . , y _(ny,rTc) _(i) )=(x _(0,rTc) _(i−1) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) , √{square root over ( α)}y _(0,rTc) _(i−1) ,y _(1,rTc) _(i−1) , . . . , x _(ny,rTc) _(i−1) ), otherwise they are computed as: (x _(0,rTc) _(i) ,x _(1,rTc) _(i) , . . . , x _(nx,rTc) _(i) ,√{square root over ( α)}y _(0,rTc) _(i) ,y _(1,rTc) _(i) , . . . , y _(ny,rTc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,√{square root over ( α)}y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rTc) _(i−1) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) ,√{square root over ( α)}y _(0,rTc) _(i−1) ,y _(1,rTc) _(i−1) , . . . , x _(ny,rTc) _(i−1) ) and (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i) , . . . , x _(ny,rc) _(i−1) ); t) the appropriate bits of the x-coordinates x_(0,rc), x_(1,rc), . . . , x_(nx,rc), y-coordinates y_(1,rc), . . . , y_(ny,rc) and the sign bit of y_(0,rc) of the message authentication code point (x_(0,rc) _(u) , x_(1,rc) _(u) , . . . , x_(nx,rc) _(u) , y_(0,rc) _(u) , y_(1,rc) _(u) , . . . , y_(ny,rc) _(u) ) and the appropriate bits of the x-coordinates x_(0,rTc), x_(1,rTc), . . . , x_(nx,rTc), y-coordinates y_(1,rTc), . . . , y_(ny,rTc) and the sign bit of y_(0,rTc) of the message authentication point (x_(0,rTc) _(u) , x_(1,rTc) _(u) , . . . , x_(nx,rTc) _(u) , √{square root over ( α)}y_(0,rTc) _(u) , y_(1,rTc) _(u) , . . . , y_(ny,rTc) _(u) ) are concatenated together to form the message authentication code of the received message; and u) if the received message authentication code is equal to the message authentication code of the received message generated at the receiving correspondent, the received message is authenticated.
 4. The computerized method of generating an elliptic polynomial-based message authentication code as recited in claim 3, wherein the non-iterative embedding method comprises the steps of: a) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), wherein M and L are integers; b) dividing each M-bit string into (nx+ny+1) strings mx₀, mx₁, . . . , mx_(nx), my₁, . . . , my_(ny), wherein the length of string mx₀ is less than or equal to (N−L) bits, the length of each of the remaining (nx+ny) strings mx₁, . . . , mx_(nx), my₁, . . . , my_(ny) being less than or equal to N bits, wherein each of the (nx+ny+1) bit strings has an equivalent value which is an element of the finite field F; c) assigning the value of the bit strings of mx₁, . . . , mx_(nx) to x₁, . . . , x_(nx); d) assigning the value of the bit strings of my₁, . . . , my_(ny) to Y₁, . . . , y_(ny); e) assigning the value of the bit string of mx₀ to x₀; f) substituting the values of x₀, x₁, . . . , x_(nx), and y₁, . . . , y_(ny) in a selected elliptic polynomial equation with (nx+1) x-coordinates and (ny+1) y-coordinates to form a quadratic equation in y₀ of the form y₀ ²+Ay₀+B=0, wherein A and B are variable coefficients; g) if the quadratic equation in y_(o) has solutions y _(0,1), y _(0,2), wherein the solutions are elements of the finite field F, then a method agreed upon by the sending and receiving correspondent is used to select one of the solutions y _(0,1), y _(0,2), and further assigning the chosen solution to y_(o), assigning α_(m)=1, and the message point is given by (x_(0,m), x_(1,m), . . . , x_(nx,m), y_(0,m), y_(1,m), . . . , y_(ny,m)), otherwise a method agreed upon by the sending and receiving correspondent is used to select one of the solutions √{square root over ( α)} y _(0,1), √{square root over ( α)} y _(0,2), and then further assign the chosen solution to y_(o), assigning α_(m)= α and the message point is given by (x_(0,m), x_(1,m), . . . , x_(nx,m), √{square root over ( α)}y_(0,m), y_(1,m), . . . , y_(ny,m)); and h) denoting the message point as (x_(0,m), x_(1,m), . . . , x_(nx,m), √{square root over (α)}y_(0,m), y_(1,m), . . . , y_(ny,m)).
 5. A computerized method of generating an elliptic polynomial-based message authentication code, comprising the steps of: a) defining a maximum block size that can be embedded into (nx+1) x-coordinates and ny y-coordinates, wherein n is an integer, and setting the maximum block size to be (nx+ny+1)N bits, wherein N is an integer; b) a sending correspondent and a receiving correspondent agree upon the values of nx and ny, and further agree on a set of coefficients a_(1k), a_(2kl), a_(3k), c_(1lki), c_(2kl), c_(3kli), b_(1l), b_(2lk), b_(3lk), b_(4k) & b_(c)εF, wherein F represents a finite field and where the elements of F are represented in N-bits, the sending and receiving correspondents further agreeing on a random number k that represents a shared secret key for communication, the sending and receiving correspondents further agreeing on a base point on an elliptic polynomial (x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B))εEC^(xn+ny+2), wherein EC^(nx+ny+2) represents the elliptic polynomial, and further agreeing on a base point on the twist of the elliptic polynomial (x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB))εEC^(xn+ny+2), wherein TEC^(nx+ny+2) represents the twist of the elliptic polynomial EC^(nx+ny+2), and wherein α is a non-quadratic residue element of F, the sending and receiving correspondents further agreeing on an initial vector s_(m) ⁻¹ having a length of (nx+ny+1)N-bits; the sending correspondent then performs the following steps: c) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny),of a key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using a non-iterative embedding method; d) if (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, a scalar multiplication is computed as (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=k(x _(0,TB) ,x _(1,TB) , . . . , x _(nx,TB) ,√{square root over ( α)}y _(0,B) ,y _(1,TB) , . . . , y _(ny,TB)) and (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,y _(0,k) ,y _(1,k) , . . . , y _(ny,k)) otherwise, if α_(k)= α, then the scalar multiplication is computed as (x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ )=k(x _(0,B) ,x _(1,B) , . . . , x _(nx,B) ,y _(0,B) ,y _(1,B) , . . . , y _(ny,B)) and (x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ )=(x _(0,k) ,x _(1,k) , . . . , x _(nx,k) ,√{square root over ( α)}y _(0,k) ,y _(1,k) , . . . , y _(ny,k)); e) denoting the first (nx+ny+1)N-bit string of the message data to be s_(m) ₀ ′ and computing the (nx+ny+1)N-bit string s_(m) ₀ as s_(m) ₀ =s_(m) ₀ ′⊕s_(m) ⁻¹ ; f) dividing the message bit string into M-bit strings, where ((nx+ny+1)N−L)>M>(N−L), wherein M and L are integers, and embedding the message bit string of the initial block of the message bit string s_(m) ₀ into the (nx+1) x-coordinates and ny y-coordinates of the elliptic message point (x_(0,m) ₀ , x_(1,m) ₀ , . . . , x_(nx,m) ₀ , √{square root over (α_(m) ₀ )}y_(0,m) ₀ , y_(1,m) ₀ , . . . , y_(ny,m) ₀ ) using the non-iterative embedding method; g) if the message point of the initial block is on the elliptic polynomial, then the message authentication code points are computed as: (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x _(0,Tc) ₀ ,x _(1,Tc) ₀ , . . . , x _(nx,Tc) ₀ , √{square root over ( α)}y _(0,Tc) ₀ ,y _(1,Tc) ₀ , . . . , y _(ny,Tc) ₀ )=(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ , √{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , y _(ny,TS) ₀ ), otherwise, they are computed as: (x _(0,Tc) ₀ ,x _(1,Tc) ₀ , . . . , x _(nx,Tc) ₀ ,√{square root over ( α)}y _(0,Tc) ₀ ,y _(1,Tc) ₀ , . . . , y _(ny,Tc) ₀ )=(x _(0,m) ₀ ,x _(1,m) ₀ , . . . , x _(nx,m) ₀ ,√{square root over ( α)}y _(0,m) ₀ ,y _(1,m) ₀ , . . . , y _(ny,m) ₀ )+(x _(0,TS) ₀ ,x _(1,TS) ₀ , . . . , x _(nx,TS) ₀ ,√{square root over ( α)}y _(0,TS) ₀ ,y _(1,TS) ₀ , . . . , x _(ny,TS) ₀ ) and (x _(0,c) ₀ ,x _(1,c) ₀ , . . . , x _(nx,c) ₀ ,y _(0,c) ₀ ,y _(1,c) ₀ , . . . , y _(ny,c) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ); h) initializing i as i=1, and repeating the following steps i) through l) and incrementing i at each step until all of the message data blocks are processed: i) denoting the (nx+ny+1)N-bit string of an l^(th) message data to be S_(m) _(i) and computing the (nx+ny+1)N-bit string s_(m) _(i) as s_(m) _(i) =s_(m) _(i) ′⊕s_(m) _(i−1) ; j) embedding the message (nx+ny+1)N-bit string of an i-th block S_(m) _(i) ′ into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx) and the ny y-coordinates, y₁, . . . , y_(ny) of the elliptic message point (x_(0,m) _(i) , x_(1,m) _(i) , . . . , x_(nx,m) _(i) , √{square root over (α_(m) _(i) )}y_(0,m) _(i) , y_(1,m) _(i) , . . . , y_(ny,m) ₀ ) using the non-iterative embedding method; k) doubling the points (x_(0,c) _(i−1) , x_(1,c) _(i−1) , . . . , x_(nx,c) _(i−1) , y_(0,c) _(i−1) , y_(1,c) _(i−1) , . . . , y_(ny,c) _(i−1) ) and (x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , y_(ny,Tc) _(i−1) ) as: (x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , y _(ny,c) _(i−1) )=2(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , y _(ny,c) _(i−1) ) and (x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , y _(ny,Tc) _(i−1) )=2(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) ,y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , y _(ny,Tc) _(i−1) ); l) if the message point of the i-th block is on the elliptic polynomial, then the message authentication code points are computed as: (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , x _(ny,c) _(i−1) ) and (x_(0,Tc) _(i) , x_(1,Tc) _(i) , . . . , x_(nx,Tc) _(i) , √{square root over ( α)}y_(0,Tc) _(i) , y_(1,Tc) _(i) , . . . , y_(ny,Tc) _(i) )=(x_(0,Tc) _(i−1) , x_(1,Tc) _(i−1) , . . . , x_(nx,Tc) _(i−1) , √{square root over ( α)}y_(0,Tc) _(i−1) , y_(1,Tc) _(i−1) , . . . , x_(ny,Tc) _(i−1) ), otherwise, they are computed as: (x _(0,Tc) _(i) ,x _(1,Tc) _(i) , . . . , x _(nx,Tc) _(i) ,√{square root over ( α)}y _(0,Tc) _(i) ,y _(1,Tc) _(i) , . . . , y _(ny,Tc) _(i) )=(x _(0,m) _(i) ,x _(1,m) _(i) , . . . , x _(nx,m) _(i) ,√{square root over ( α)}y _(0,m) _(i) ,y _(1,m) _(i) , . . . , y _(ny,m) _(i) )+(x _(0,Tc) _(i−1) ,x _(1,Tc) _(i−1) , . . . , x _(nx,Tc) _(i−1) , √{square root over ( α)}y _(0,Tc) _(i−1) ,y _(1,Tc) _(i−1) , . . . , x _(ny,Tc) _(i−1) ) and (x _(0,c) _(i) ,x _(1,c) _(i) , . . . , x _(nx,c) _(i) ,y _(0,c) _(i) ,y _(1,c) _(i) , . . . , y _(ny,c) _(i) )=(x _(0,c) _(i−1) ,x _(1,c) _(i−1) , . . . , x _(nx,c) _(i−1) ,y _(0,c) _(i−1) ,y _(1,c) _(i−1) , . . . , x _(ny,c) _(i−1) ); m) a set of appropriate bits of the x-coordinates x_(0,c), x_(1,c), . . . , x_(nx,c), y-coordinates y_(1,c), . . . , y_(ny,c) and a sign bit of y_(0,c) of the message authentication code point (x_(0,c) _(u) , x_(1,c) _(u) , . . . , x_(nx,c) _(u) , y_(0,c) _(u) , y_(1,c) _(u) , . . . , y_(ny,c) ₀ ), and a set of appropriate bits of the x-coordinates x_(0,Tc), x_(1,Tc), . . . , x_(nx,Tc), y-coordinates y_(1,Tc), . . . , y_(ny,Tc) and the sign bit of y_(0,Tc) of the message authentication code point (x_(0,Tc) _(u) , x_(1,Tc) _(u) , . . . , x_(nx,Tc) _(u) , √{square root over ( α)}y_(0,Tc) _(u) , y_(1,Tc) _(u) , . . . , y_(ny,Tc) _(u) ) are concatenated together to form the message authentication code, which is then appended to the message to be sent together with, if necessary, any additional information needed to help authenticate the message bit string at the receiving correspondent without compromising security and sent to the receiving correspondent; the receiving correspondent then performs the following steps: n) embedding the bit string of the secret key into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of the key elliptic point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) using the non-iterative embedding method; o) if the point (x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over (α_(k))}y_(0,k), y_(1,k), . . . , y_(ny,k)) is on the elliptic polynomial, then a scalar multiplication is performed as (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=k(x_(0,TB), x_(1,TB), . . . , x_(nx,TB), √{square root over ( α)}y_(0,B), y_(1,TB), . . . , y_(ny,TB)) and (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), y_(0,k), y_(1,k), . . . , y_(ny,k)), otherwise if α_(k)=α_(o) then the scalar multiplication is performed as (x_(0,S) ₀ , x_(1,S) ₀ , . . . , x_(nx,S) ₀ , y_(0,S) ₀ , y_(1,S) ₀ , . . . , x_(ny,S) ₀ )=k(x_(0,B), x_(1,B), . . . , x_(nx,B), y_(0,B), y_(1,B), . . . , y_(ny,B)) (x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ )=(x_(0,k), x_(1,k), . . . , x_(nx,k), √{square root over ( α)}y_(0,k), y_(1,k), . . . , y_(ny,k)); p) denoting the first (nx+ny+1)N-bit string of the message data to be s_(rm) ₀ ′ and computing the (nx+ny+1)N-bit string s_(rm) ₀ as s_(rm) ₀ n=s_(rm) ₀ ′⊕s_(rm) ⁻¹ ; q) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), wherein M and L are integers, and embedding the message bit string of the initial block of the message bit string s_(rm) ₀ into the (nx+1) x-coordinates and the ny y-coordinates of the received message elliptic point (x_(0,rm) ₀ , x_(1,rm) ₀ , . . . , x_(nx,rm) ₀ , √{square root over (α_(rm) ₀ )}y_(0,rm) ₀ , y_(1,rm) ₀ , . . . , y_(ny,rm) ₀ ) using the non-iterative embedding method; r) if the message point of the 0-th block is on the elliptic-polynomial, then the message authentication points are computed as: (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ) and (x_(0,rTc) ₀ , x_(1,rTc) ₀ , . . . , x_(nx,rTc) ₀ , √{square root over ( α)}y_(0,rTc) ₀ , y_(1,rTc) ₀ , . . . , y_(ny,rTc) ₀ )=(x_(0,TS) ₀ , x_(1,TS) ₀ , . . . , x_(nx,TS) ₀ , √{square root over ( α)}y_(0,TS) ₀ , y_(1,TS) ₀ , . . . , x_(ny,TS) ₀ ) otherwise they are computed as: (x _(0,rTc) ₀ ,x _(1,rTc) ₀ , . . . , x _(nx,rTc) ₀ ,√{square root over ( α)}y _(0,rTc) ₀ ,y _(1,rTc) ₀ , . . . , y _(ny,rTc) ₀ )=(x _(0,rm) ₀ ,x _(1,rm) ₀ , . . . , x _(nx,rm) ₀ ,√{square root over ( α)}y _(0,rm) ₀ ,y _(1,rm) ₀ , . . . , y _(ny,rm) ₀ )+(x _(0,rTS) ₀ ,x _(1,rTS) ₀ , . . . , x _(nx,rTS) ₀ ,√{square root over ( α)}y _(0,rTS) ₀ ,y _(1,rTS) ₀ , . . . , x _(ny,rTS) ₀ ) and (x _(0,rc) ₀ ,x _(1,rc) ₀ , . . . , x _(nx,rc) ₀ ,y _(0,rc) ₀ ,y _(1,rc) ₀ , . . . , y _(ny,rc) ₀ )=(x _(0,S) ₀ ,x _(1,S) ₀ , . . . , x _(nx,S) ₀ ,y _(0,S) ₀ ,y _(1,S) ₀ , . . . , x _(ny,S) ₀ ); s) initializing the integer i as i=0, and repeating the following steps t) through w) and incrementing i at each step until all of the message data blocks are processed: t) denoting the (nx+ny+1)N-bit string of an l^(th) message data to be s_(rm) _(i) ′ and computing the (nx+ny+1)N-bit string s_(rm) _(i) as s_(rm) _(i) =s_(rm) _(i) ′⊕s_(rm) _(i−1) ; u) embedding the received message (nx+ny+1)N-bit string of the i-th block s_(rm) _(i) the into the (nx+1) x-coordinates x₀, x₁, . . . , x_(nx), and the ny y-coordinates y₁, . . . , y_(ny), of the received message elliptic point (x_(0,rm) _(i) , x_(1,rm) _(i) , . . . , x_(nx,rm) _(i) , √{square root over (α_(rm) _(i) )}y_(0,rm) _(i) , y_(1,rm) _(i) , . . . , y_(ny,rm) _(i) ) using the non-iterative embedding method; v) doubling the points (x_(0,rc) _(i−1) , x_(1,rc) _(i−1) , . . . , x_(nx,rc) _(i−1) , y_(0,rc) _(i−1) , y_(1,rc) _(i−1) , . . . , y_(ny,rc) _(i−1) ) and (x_(0,rTc) _(i−1) , x_(1,rTc) _(i−1) , . . . , x_(nx,rTc) _(i−1) , y_(0,rTc) _(i−1) , y_(1,rTc) _(i−1) , . . . , y_(ny,rTc) _(i−1) ) as: (x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i−1) , . . . , y _(ny,rc) _(i−1) )=2(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i−1) , . . . , y _(ny,rc) _(i−1) ) (x _(0,rTc) _(i−1) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) ,y _(0,rTc) _(i−1) ,y _(1,rTc) _(i−1) , . . . , y _(ny,rTc) _(i−1) )=2(x _(0,rTc) _(i−1) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) ,y _(0,rTc) _(i−1) ,y _(1,rTc) _(i−1) , . . . , y _(ny,rTc) _(i−1) ); w) if the received message point of the i-th block is on the elliptic polynomial, then the message authentication points are computed as (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i) , . . . , x _(ny,rc) _(i−1) ) and (x_(0,rTc) _(i) , x_(1,rTc) _(i) , . . . , x_(nx,rTc) _(i) , √{square root over ( α)}y_(0,rTc) _(i) , y_(1,rTc) _(i) , . . . , y_(ny,rTc) _(i) )=(x_(0,rTc) _(i−1) , x_(1,rTc) _(i−1) , . . . , x_(nx,rTc) _(i−1) , √{square root over ( α)}y_(0,rTc) _(i−1) , y_(1,rTc) _(i−1) , . . . , x_(ny,rTc) _(i−1) ), otherwise, they are computed as: (x _(0,rTc) _(i) ,x _(1,rTc) _(i) , . . . , x _(nx,rTc) _(i) ,√{square root over ( α)}y _(0,rTc) _(i) ,y _(1,rTc) _(i) , . . . , y _(ny,rTc) _(i) )=(x _(0,rm) _(i) ,x _(1,rm) _(i) , . . . , x _(nx,rm) _(i) ,√{square root over ( α)}y _(0,rm) _(i) ,y _(1,rm) _(i) , . . . , y _(ny,rm) _(i) )+(x _(0,rTc) _(i−1) ,x _(1,rTc) _(i−1) , . . . , x _(nx,rTc) _(i−1) ,√{square root over ( α)}y _(0,rTc) _(i−1) ,y _(1,rTc) _(i−1) , . . . , x _(ny,rTc) _(i−1) ) and (x _(0,rc) _(i) ,x _(1,rc) _(i) , . . . , x _(nx,rc) _(i) ,y _(0,rc) _(i) ,y _(1,rc) _(i) , . . . , y _(ny,rc) _(i) )=(x _(0,rc) _(i−1) ,x _(1,rc) _(i−1) , . . . , x _(nx,rc) _(i−1) ,y _(0,rc) _(i−1) ,y _(1,rc) _(i) , . . . , x _(ny,rc) _(i−1) ); x) the appropriate bits of the x-coordinates x_(0,rc), x_(1,rc), . . . , x_(nx,rc), y-coordinates y_(1,rc), . . . , y_(ny,rc) and the sign bit of y_(0,rc) of the message authentication code point (x_(0,rc) _(u) , x_(1,rc) _(u) , . . . , x_(nx,rc) _(u) , y_(0,rc) _(u) , y_(1,rc) _(u) , . . . , y_(ny,rc) ₀ ) and the appropriate bits of the x-coordinates x_(0,rTc), x_(1,rTc), . . . , x_(nx,rTc), y-coordinates y_(1,rTc), . . . , y_(ny,rTc) and the sign bit of y_(0,rTc) of the message authentication point (x_(0,rTc) _(u) , x_(1,rTc) _(u) , . . . , x_(nx,rTc) _(u) , √{square root over ( α)}y_(0,rTc) _(u) , y_(1,rTc) _(u) , . . . , y_(ny,rTc) _(u) ) are concatenated together to form the message authentication code of the received message; and y) if the received message authentication code is equal to the message authentication code of the received message generated at the receiving correspondent, the received message is authenticated.
 6. The computerized method of generating an elliptic polynomial-based message authentication code as recited in claim 5, wherein the non-iterative embedding method comprises the steps of: a) dividing the message bit string into M-bit strings where ((nx+ny+1)N−L)>M>(N−L), wherein M and L are integers; b) dividing each M-bit string into (nx+ny+1) strings mx₀, mx₁, . . . , mx_(nx), my₁, . . . , my_(ny), wherein the length of string mx₀ is less than or equal to (N−L) bits, the length of each of the remaining (nx+ny) strings mx₁, . . . , mx_(nx), my₁, . . . , my_(ny) being less than or equal to N bits, wherein each of the (nx+ny+1) bit strings has an equivalent value which is an element of the finite field F; c) assigning the value of the bit strings of mx₁, . . . , mx_(nx) to x₁, . . . , x_(nx); d) assigning the value of the bit strings of my₁, . . . , my_(ny) to y₁ . . . y_(ny); e) assigning the value of the bit string of mx₀ to x₀; f) substituting the values of x₀, x₁, . . . , x_(nx) and y₁, . . . , y_(ny) in a selected elliptic polynomial equation with (nx+1) x-coordinates and (ny+1) y-coordinates to form a quadratic equation in y₀ of the form y₀ ²+Ay₀+B=0, wherein A and B are variable coefficients; g) if the quadratic equation in y_(o) has solutions y _(0,1), y _(0,2), wherein the solutions are elements of the finite field F, then a method agreed upon by the sending and receiving correspondent is used to select one of the solutions y _(0,1), y _(0,2), and further assigning the chosen solution to y_(o), assigning α_(m)=1, and the message point is given by (x_(0,m), x_(1,m), . . . , x_(nx,m), y_(0,m), y_(1,m), . . . , y_(ny,m)), otherwise a method agreed upon by the sending and receiving correspondent is used to select one of the solutions √{square root over ( α)} y _(0,1), √{square root over ( α)} y _(0,2), and then further assign the chosen solution to y_(o), assigning α_(m)= α and the message point is given by (x_(0,m), x_(1,m), . . . , x_(nx,m), √{square root over ( α)}y₀, y_(1,m), . . . , y_(ny,m)); and h) denoting the message point as (x_(0,m), x_(1,m), . . . , x_(nx,m), √{square root over (α)}y_(0,m), y_(1,m), . . . , y_(ny,m)). 